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5.1 The Concept of Decoherence 

This introduction to the theory of decoherence is aimed at readers with an 
interest in the science of quantum information. In that field, one is usually 
content with simple, abstract descriptions of non-unitary "quantum chan- 
nels" to account for imperfections in quantum processing tasks. However, in 
order to justify such models of non-unitary evolution and to understand their 
limits of applicability it is important to know their physical basis. I will there- 
fore emphasize the dynamic and microscopic origins of the phenomenon of 
decoherence, and will relate it to concepts from quantum information where 
applicable, in particular to the theory of quantum measurement. Q 

The study of decoherence, though based at the heart of quantum theory, 
is a relatively young subject. It was initiated in the 1970's and 1980's with 
the work of H. D. Zeh and W. Zurek on the emergence of classicality in the 
quantum framework. Until that time the orthodox interpretation of quantum 
mechanics dominated, with its strict distinction between the classical macro- 
scopic world and the microscopic quantum realm. The mainstream attitude 
concerning the boundary between the quantum and the classical was that 
this was a purely philosophical problem, intangible by any physical analysis. 
This changed with the understanding that there is no need for denying quan- 
tum mechanics to hold even macroscopically, if one is only able to understand 
within the framework of quantum mechanics why the macro-world appears 
to be classical. For instance, macroscopic objects are found in approximate 
position eigenstates of their center-of-mass, but never in superpositions of 
macroscopically distinct positions. The original motivation for the study of 
decoherence was to explain these effective super- selection rules and the ap- 
parent emergence of classicality within quantum theory by appreciating the 
crucial role played by the environment of a quantum system. 

Hence, the relevant theoretical framework for the study of decoherence 
is the theory of open quantum systems, which treats the effects of an un- 
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controllable environment on the quantum evolution. Originally developed to 
incorporate the phenomena of friction and thermalization in the quantum 
formulation, it has of course a much longer history than decoherence theory. 
However, we will see that the intuition and approximations developed in the 
traditional treatments of open quantum systems are not necessarily appropri- 
ate to yield a correct description of decoherence effects, which may take place 
on a time scale much shorter than typical relaxation phenomena. In a sense, 
one may say that while the traditional treatments of open quantum systems 
focus on how an environmental "bath" affects the system, the emphasis in 
decoherence is more on the contrary question, namely how the system affects 
and disturbs environmental degrees of freedom, thereby revealing information 
about its state. 

The physics of decoherence became very popular in the last decade, mainly 
due to advances in experimental technology. In a number of experiments 
the gradual emergence of classical properties in a quantum system could be 
observed, in agreement with the predictions of decoherence theory. Needles 
to say, a second important reason for the popularity of decoherence is its 
relevance for quantum information processing tasks, where the coherence of 
a relatively large quantum system has to be maintained over a long time. 

Parts of these lecture notes are based on the books on decoherence by 
E. Joos et al. [T] and on open quantum systems by H.-P. Breuer & F. Petruc- 
cione [2], and on the lecture notes of W. Strunz [3]. Interpretational aspects, 
which are not covered here, are discussed in [H [5], and useful reviews by 
W. H. Zurek and J. P. Paz can be found in [6l[7] . This chapter deals exclusively 
with conventional, i.e. environmental decoherence, as opposed to spontaneous 
reduction theories [8], which aim at "solving the measurement problem of 
quantum mechanics" by modifying the Schrodinger equation. These models 
are conceptually very different from environmental decoherence, though their 
predictions of super-selection rules are often qualitatively similar. 

5.1.1 Decoherence in a Nutshell 

Let us start by discussing the basic decoherence effect in a rather general 
framework. As just mentioned, we need to account for the unavoidable cou- 
pling of the quantum system to its environment. Although these environmen- 
tal degrees of freedom are to be treated quantum mechanically, their state 
must be taken unobservable for all practical purposes, be it due to their large 
number or uncontrollable nature. In general, the detailed temporal dynam- 
ics induced by the environmental interaction will be very complicated, but 
one can get an idea of the basic effects by assuming that the interaction is 
sufficiently short-ranged to admit a description in terms of scattering the- 
ory. In this case, only the map between the asymptotically free states before 
and after the interaction needs to be discussed, thus avoiding a temporal 
description of the collision dynamics. 
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Let the quantum state of a system be described by the density opera- 
tor p on the Hilbert space Tt. We take the system to interact with a single 
environmental degree of freedom at a time - think of a phonon, a polaron, 
or a gas particle scattering off your favorite implementation of a quantum 
register. Moreover, let us assume, for the time being, that this environmental 
"particle" is in a pure state pe = |'0in)(V'in|E, with \ipi n )E £ He- The scat- 
tering operator Stot maps between the in- and out-asymptotes in the total 
Hilbert space Wtot = H <£> Hb, and for sufficiently short-ranged interaction 
potentials we may identify those with the states before and after the collision. 
The initially uncorrelated system and environment turn into a joint state, 

[before collision] p to t = p <%> \ipia) , (5.1) 

[after collision] p' tot = S tot [p® |V'in)(V'in|E]Sj ot . (5.2) 

Now let us assume, in addition, that the interaction is non-invasive with 
respect to a certain system property. This means that there is a number 
of distinct system states, such that the environmental scattering off these 
states causes no transitions in the system. For instance, if these distinguished 
states correspond to the system being localized at well-defined sites then the 
environmental particle should induce no hopping between the locations. In 
the case of elastic scattering, on the other hand, they will be the energy 
eigenstates. Denoting the set of these mutually orthogonal system states by 
{\n)} £ Ti, the requirement of non-invasiveness means that Stot commutes 
with those states, that is, it has the form 

Stot = ^2 \ n )( n \ ® s « ' ( 5 - 3 ) 

n 

where the S„ are scattering operators acting in the environmental Hilbert 
space. The insertion into (|5.2|) yields 

Ptot = ~J2( m \p\n)\m)(n\ <g> S m |^in)(V>in|ESj, 

m,n 

and disregarding the environmental state by performing a partial trace we 
get the system state after the interaction: 

p = tr E (Pto t ) = |m)(n|p m „(^in| S^S m |Vin)E • (5.5) 

S v ' 

Since the S n are unitary the diagonal elements, or populations, are indeed 
unaffected, 



Pram ~ Pram , 



(5.6) 
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while the off-diagonal elements, or coherences, get multiplied by the overlap 
of the environmental states scattered off the system states m and n, 

P' m n=Pmn{^ t \^)- (5-7) 

This factor has a modulus of less than one so that the coherences, which 
characterize the ability of the system state to display a superposition be- 
tween |m) and get suppressed^ It is important to note that this loss of 
coherence occurs in a special basis, which is determined only by the scattering 
operator, i.e. by the type of environmental interaction, and to a degree that 
is determined by both the environmental state and the interaction. 

This loss of the ability to show quantum behavior due to the interaction 
with an environmental quantum degree of freedom is the basic effect of deco- 
herence. One may view it as due to the arising correlation between the system 
with the environment. After the interaction the joint quantum state of system 
and environment is no longer separable, and part of the coherence initially 
located in the system now resides in the non-local correlation between system 
and the environmental particle; it is lost once the environment is disregarded. 
A complementary point of view argues that the interaction constitutes an in- 
formation transfer from the system to the environment. The more the overlap 
in (|5.7[) differs in magnitude from unity, the more an observer could in prin- 
ciple learn about the system state by measuring the environmental particle. 
Even though this measurement is never made, the complementarity princi- 
ple then explains that the wave-like interference phenomenon characterized 
by the coherences vanishes as more information discriminating the distinct, 
"particle-like" system states is revealed. 

To finish the introduction, here is a collection of other characteristics and 
popular statements about the dccohcrcncc phenomenon. One often hears that 
decoherence (i) can be extremely fast as compared to all other relevant time 
scales, (ii) that it can be interpreted as an indirect measurement process, a 
monitoring of the system by the environment, (Hi) that it creates dynam- 
ically a set of preferred states ("robust states" or "pointer states") which 
seemingly do not obey the superposition principle, thus providing the mi- 
croscopic mechanism for the emergence of effective super-selection rules, and 
(iv) that it helps to understand the emergence of classicality in a quantum 
framework. These points will be illustrated in the following, though (Hi) has 
been demonstrated only for very simple model systems and (iv) depends to 
a fair extent on your favored interpretation of quantum mechanics. 

5.1.2 General Scattering Interaction 

In the above demonstration of the decoherence effect the choice of the in- 
teraction and the environmental state was rather special. Let us therefore 

1 The value \p mn \ determines the maximal fringe visibility in a general interference 
experiment involving the states \m) and |n), as described by the projection on a 
general superposition \ipe,ip) = cos (9) \m)+e' lp sin(9) \n). 
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now take S to t aud Pe to be arbitrary and carry out the same analysis. Per- 
forming the trace in (|5-5[) in the eigenbasis of the environmental state, 
Pe = J2e PM>t) {i>t\ e, we have 

p' = tr E (S to t[p ® PE]Stot) = W(V'j|Stot|'0f)EP(V'£|St t otl^)E 

= ^W feP Wt, (5.8) 

k 

where the (ipj\Stot\tpe) are operators in TL. After subsuming the two indices 
j, £ into a single one, we get the second line with the Kraus operators 
given by 

W fe = V%(ViJS tot |^J . (5.9) 
It follows from the unitarity of Stot that they satisfy 

5^w£w fc =I. (5.10) 

fc 

This implies that (|5.8[) is the operator-sum representation of a completely 
positive map <P : p i— > p' (see Sect. I5.3.1[) . In other words, the scattering 
transformation has the form of the most general evolution of a quantum 
state that is compatible with the rules of quantum theory. In the operational 
formulation of quantum mechanics this transformation is usually called a 
quantum operation [9], the quantum information community likes to call it 
a quantum channel. Conversely, given an arbitrary quantum channel, one 
can also construct a scattering operator Stot and an environmental state pe 
giving rise to the transformation, though it is usually not very helpful from a 
physical point of view to picture the action of a general, dissipative quantum 
channel as due to a single scattering event. 

5.1.3 Decoherence as an Environmental Monitoring Process 

We are now in a position to relate the decoherence of a quantum system to 
the information it reveals to the environment. Since the formulation is based 
on the notion of an indirect measurement it is necessary to first collect some 
aspects of measurement theory [10l E] . 

Elements of General Measurement Theory 

Projective Measurements 

This is the type of measurement discussed in standard textbooks of quantum 
mechanics. A projective operator |a)(a| = P a = P 2 a = is attributed to 
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each possible outcome a of an idealized measurement apparatus. The prob- 
ability of the outcome a is obtained by the Born rule 

Prob(a|p) = tr (P a p) = (a\p\a) , (5.11) 

and after the measurement of a the state of the quantum system is given by 
the normalized projection 

M --^ M ^=^S)- (5 - 12) 

The basic requirement that the projectors form a resolution of the identity 
operator, 

5^P a =I, (5.13) 

a 

ensures the normalization of the corresponding probability distribution 
Prob (a\p). 

If the measured system property corresponds to a self-adjoint operator A 
the P Q are the projectors into its eigenspaces, so that its expectation value is 

(A) = tr (Ap) . 

If A has a continuous spectrum the outcomes are characterized by intervals 
of a real parameter, and the sum in (|5.13|) should be replaced by a projector- 
valued Stieltjes integral J dP(a) =1, or equivalently by a Lebesgue integral 
over a projector-valued measure (PVM) [101 111] . 

It is important to note that projective measurements are not the most 
general type of measurement compatible with the rules of quantum mechan- 
ics. In fact, non-destructive measurements of a quantum system are usually 
not of the projective kind. 

Generalized Measurements 

In the most general measurement situation, a positive (and therefore hermi- 
tian) operator F Q > is attributed to each outcome a. Again, the collection 
of operators corresponding to all possible outcomes must form a resolution 
of the identity operator, 

5^F a =1. (5.14) 

a 

In particular, one speaks of a positive operator-valued measure (POVM) in 
the case of a continous outcome parameter, J dF(a) = I , and the probability 
(or probability density in the continuous case) of outcome a is given by 



Prob(a|p) = tr(F a p) 



(5.15) 
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The effect on the system state of a generalized measurement is described by 
a nonlinear transformation 

Al : p ^ Af (p Q ) = ' (5.16) 

tr (F a p) 

involving a norm-decreasing completely positive map in the numerator (see 
Sect. [STBTTj) , and a normalization which is subject to the consistency require- 
ment 



E 1 



M Q , fe = F a . (5.17) 



The operators M Qj fc appearing in (I5.16|) are called measurement operators, 
and they serve to characterize the measurement process completely. The F a 
are sometimes called "effects" or "measurement elements" . Note that different 
measurement operators M a ^ can lead to the same measurement element F Q . 

A simple class of generalized measurements are unsharp measurements, 
where a number of projective operators are lumped together with proba- 
bilistic weights in order to account for the finite resolution of a measurement 
device or for classical noise in its signal processing. However, generalized mea- 
surements schemes may also perform tasks which are seemingly impossible 
with a projective measurement, such as the error-free discrimination of two 
non-orthogonal states p~21 [13] . 

Efficient Measurements 

A generalized measurement is called efficient if there is only a single summand 
in (|5. 16|) for each outcome a, 

, (6 . 18) 

tr [UlU aP j 

implying that pure states are mapped to pure states. In a sense, these are 
measurements where no unnecessary, that is no classical, uncertainty is intro- 
duced during the measurement process, see below. By means of a (left) polar 
decomposition and the consistency requirement (|5. 17|) efficient measurement 
operators have the form 

M Q = U Q yF^, (5.19) 

with an unitary operator U a . This way the state after efficient measurement 
can be expressed in a form which decomposes the transformation into a "raw 
measurement" described by the F Q and a "measurement back-action" given 
by the U Q : 

M(p\a)= ^(F^pf ^ ■ (5 - 20) 

back-action * v ' back-action 

raw 
measurement 
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In this transformation the positive operators y/WZ "squeeze" the state 
along the measured property and expand it along the other, complementary 
ones, similar to what a projector would do, while the back-action operators U a 
"kick" the state by transforming it in a way that is reversible, in principle, 
provided the outcome a is known. Note that the projective measurements 
(|5.12[) are a sub-class in the set of back-action-free efficient measurements. 



Indirect Measurements 



In an indirect measurement one tries to obtain information about the system 
in a way that disturbs it as little as possible. This is done by letting a well- 
prepared microscopic quantum probe interact with the system. This probe is 
then measured by projection, i.e. destructively, so that one can infer proper- 
ties of the system without having it brought into contact with a macroscopic 
measurement device. Let p plo be be the prepared state of the probe, Stot de- 
scribe the interaction between system and probe, and P a be the projectors 
corresponding to the various outcomes of the probe measurement. The prob- 
ability of measuring a is determined by the reduced state of the probe after 
interaction, i.e. 

Prob(a|p) = tr prob c (PaP pr obo) = tr pro bc (P a tr sys (S tot [p (g> p pr obe]Sj ot )J . 

(5.21) 

By pulling out the system trace (extending the projectors to Htot = "H (8 H p ) 
and using the cyclic permutability of operators under the trace we have 

ProbHp) = tr (s\ ot [I $ P Q ]S to t [p <8 /Vobe]) = tr (F a p) , (5.22) 

with microscopically defined measurement elements 

F Q = tr probc (s t f ot [I (g> PJStot [I O Pprobc]) > (5.23) 

satisfying J2 a F Q = I . Since the probe measurement is projective, we can 
also specify the new system state conditioned on the click at a of the probe 
detector, 

M (p\a) = tr probc (Mtot (ptot|a)) 

[I g Pq]Stot [P ® Pprobc] S| ot [I Pg] \ 



tr 



probe 



tr(F Q p) 



= V^wJf- (5 ' 24) 

In the last step a convex decomposition of the initial probe state into pure 
states was inserted, p pro bc = J2k w k\' t l J k)(4>k\ ■ Taking P a — |a)(a| we thus 
get a microscopic description also of the measurement operators, 
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M 



ot,k 



y/w^(a\S to t\tpk} • 



(5.25) 



This shows that an indirect measurement is efficient (as defined above) if the 
probe is initially in a pure state, i.e. if there is no uncertainty introduced 
in the measurement process, apart from the one imposed by the uncertainty 
relations on p pro bc- 

If we know that an indirect measurement has taken place, but do not know 
its outcome a we have to resort to a probabilistic (Bayesian) description of the 
new system state. It is given by the sum over all possible outcomes weighted 
by their respective probabilities, 



This form is the same as above in (|5.8|) and (15. 9j) . where the basic effect of 
decoherence has been described. This indicates that the decoherence process 
can be legitimately viewed as a consequence of the information transfer from 
the system to the environment. The complementarity principle can then be 
invoked to understand which particular system properties lose their quantum 
behavior, namely those complementary to the ones revealed to the environ- 
ment. This "monitoring interpretation" of the decoherence process will help 
us below to derive microscopic master equations. 

5.1.4 A Few Words on Nomenclature 

Since decoherence phenomena show up in quite different sub-communities of 
physics, a certain confusion and lack of uniformity developed in the termi- 
nology. This is exacerbated by the fact that decoherence often reveals itself 
as a loss of fringe visibility in interference experiments - a phenomenon, 
though, which may have other causes than decoherence proper. Here is an 
attempt of clarification: 

decoherence: In the original sense, an environmental quantum effect af- 
fecting macroscopically distinct states. The term is nowadays applied to 
mesoscopically different states as well, and even for microscopic states, as 
long as it refers to the quantum effect of environmental, i.e. in practice 
unobservable, degrees of freedom. 

However, the term is often (ab-)used for any other process reducing the 
purity of a micro-state. 

dephasing: In a narrow sense, this describes the phenomenon that coher- 
ences, i.e., the off-diagonal elements of the density matrix, get reduced 
in a particular basis, namely the energy eigenbasis of the system. It is 
a statement about the effect, and not the cause. In particular, dephas- 
ing may be reversible if it is not due to decoherence, as revealed e.g. in 
spin-echo experiments. 




(5.26) 
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This phrase should be treated with great care since it is used differently 
in various sub-communities. It is taken as a synonym to " dispersion''' in 
molecular physics and in nonlinear optics, as a synonym to " decoherence" 
in condensed matter, and often as a synonym to "phase averaging" in 
matter wave optics. It is also called a T^-process in NMR and in condensed 
matter physics (see below). 

phase averaging: A classical noise phenomenon entering through the de- 
pendence of the unitary system evolution on external control parameters 
which fluctuate (parametric average over unitary evolutions). A typical 
example are the vibrations of an interferometer grating or the fluctua- 
tions of the classical magnetic field in an electron interferometer due to 
technical noise. Empirically, phase averaging is often hard to distinguish 
from decoherence proper. 

dispersion: Coherent broadening of wave packets during the unitary evo- 
lution, e.g. due to a velocity dependent group velocity or non-harmonic 
energy spacings. This unitary effect may lead to a reduction of signal 
oscillations, for instance, in molecular pump-probe experiments. 
dissipation: Energy exchange with the environment leading to thermal- 
ization. Usually accompanied by decoherence, but see Sect. 15.3.41 for a 
counter-example. 

5.2 Case Study: Dephasing of Qubits 

So far, the discussion of the temporal dynamics of the decoherence process 
was circumvented by using a scattering description. Before going to the gen- 
eral treatment of open quantum systems in Sect. 15.31 it is helpful to take a 
closer look on the time evolution of a special system where the interaction 
with a model environment can be treated exactly [HJ [2] . 

5.2.1 An Exactly Solvable Model 

Let us take a two-level system, or qubit, described by the Pauli spin operator 
a z , and model the environment as a collection of bosonic field modes. In 
practice, such fields can yield an appropriate effective description even if the 
actual environment looks quite differently in particular if the environmental 
coupling is a sum of many small contributionso What is fairly non-generic in 
the present model is the type of coupling between system and environment, 
which is taken to commute with the system Hamiltonian. 
The total Hamiltonian thus reads 



2 A counter-example would be the presence of a degenerate environmental degree 
of freedom, such as a bistable fluctuator. 
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Htot = -=-o-, + ^ b I b fe +^ (s* b l + 9k bk ) ' ( 5 - 27 ) 



Ho Hi n t 

with the usual commutation relation for the mode operators of the bosonic 
field modes, [bj, bt] = So-, and coupling constants gk- The fact that the 
system Hamiltonian commutes with the interaction, guarantees that there is 
no energy exchange between system and environment so that we expect pure 
dephasing. 

By going into the interaction picture one transfers the trivial time evo- 
lution generated by Ho to the operators (and indicates this with a tilde). In 
particular, 

Hi„t (t) = e^^Hinte-^'^ = o z ^(^e^b* + ^e^b*) , (5.28) 

k 

where the second equality is granted by the commutation [a z , H; nt ] = 0. The 
time evolution due to this Hamiltonian can be formally expressed as a Dyson 
series, 

U(t) = T_ exp (~ J dt'H mt (t') ) 

00 1 / 1 \ n /"* r~ ~ 1 

= 22—[ r n) / d * 1 "' dtnT <- [ H int(<l)---H in t(irOj . (5-29) 
n=0 " V 1 / J® 

where is the time ordering operator (putting the operators with larger 
time arguments to the left) . Due to this time ordering requirement the series 
usually cannot be evaluated exactly (if it converges at all). However, in the 
present case the commutator of Hi nt at different times is not an operator, but 
just a c-number, 

[Hi„t W,H int (i')] = \9k\ 2 sin ( Lo k (t' -t)) . (5.30) 

k 

As a consequence, the time evolution differs only by a time-dependent phase 
from the one obtained by casting the operators in their natural order|f| 

U(t) = e^W exp (-i jf dt'H int (i')) , (5.31) 



3 To obtain the time evolution U (t) for the case [H (t) , H (t')] — cl define the 
operators 



*(t) = i| < dt'H(t') 



and U (t) — exp[i<t> (t)]U (t). This way U (t) describes the "additional" motion 
due to the time ordering requirement. It satisfies 
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where the phase is given by 



Hint (tl) , Hint (£2) 



(5.32) 



^W = ^| dti jf dt 2 0(h-t 2 ) 
One can now perform the integral over the interaction Hamiltonian to get 

U(t) - e^» exp (1^ Yl ( a k(t)b{ a* k (t)b k ) ) , (5.33) 



with complex, time-dependent functions 



a k (t) := 2g k 



frw k 



(5.34) 



The operator U(i) is diagonal in the eigenbasis of the system, and it describes 
how the environmental dynamics depends on the state of the system. In 
particular, if the system is initially in the upper level, = | f), one has 



U(t)| T>|€o>e = e 1 ^! T) II Da 



a k (t) 



l£o> =:e^W| T)|^ T (*)>e , (5.35) 



and for the lower state 



U(i)|l)|&)E = ^ , »|l)nD»(-^f)|6) =: e^lD^^E. 

(5.36) 

Here we introduced the unitary displacement operators for the fc-th field 
mode. 



D k (a) = exp(ab{ - a*b k ) , 



(5.37) 



dtU(t) 



d_ i4 »(t) 
dt 



e -i* (t)+ li»(t)fj (t)e -i»W 



U(t) 



The derivative in square brackets has to be evaluated with care since the H (t) do 
not commute at different times. By first showing that [A, <9tA] = c G C implies 
d t A n = nA n - 1 a t A - \n (n - 1) cA n ~ 2 one finds 



-1 i*(t)p 



t2\-i r« 



i*(t) 



dt'H (*') , H (t) 



Therefore, we have dtU (t) = (2ft 2 ) f Q dt'[H (*') , H (t)] U (t), which can be in- 
tegrated to yield finally 



U (t) = exp JL jf* dt! jf** dfa [H («0 , H (fa)] 



i«(t) 
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which effect a translation of the field state in its attributed phase space. In 
particular, the coherent state \a) k of the field mode k is obtained from its 
ground state \0) k by \a) k := D fc (a) |0) fc [15]. 

The equations (|5.35|) and (|5.36p show that the collective state of the field 
modes gets displaced by the interaction with the system and that the sense 
of the displacement is determined by the system state. 

Assuming that the states of system and environment are initially uncor- 
rected, ptot (0) = p® p E , the time-evolved system state reads0 

p(t)=tr E ([J(t)[ P ®p E p(t^ . (5.38) 

It follows from (|5 .35[) and (|5.36p that the populations are unaffected, 

(T \Kt)\ T> = (T l/5(0)| t) , 
U \p(t)\ i) = U Ip(0)| 4) , 

while the coherences are suppressed by a factor which is given by the trace 
over the displaced initial field state, 

<T 4) = (T l/5(0)| |) tr E (J] D k (a k (t)) PE j . (5.39) 

" v ' 

x(t) 

Incidentally, the complex suppression factor \ (t) is equal to the Wigner char- 
acteristic function of the original environmental state at the points a k (t), i.e. 
it is given by the Fourier transform of its Wigner function [16] . 



Initial Vacuum State 



If the environment is initially in its vacuum state, p E = (& k |0)(0|fc, the 
|£j (t)} E and |£j.(i))E defined in (I5.35p . (I5.36P turn into multi-mode coherent 
states, and the suppression factor can be calculated immediately to yield: 

XvacW = H(0\D k (a k (t))\0) k = J^cxp ( -MC j 

k k V / 

-x P f-E4| 3fe | 2 ^S^l ■ (5-40) 



For times that are short compared to the field dynamics, t <C , one 
observes a Gaussian decay of the coherences. Modifications to this become 



4 In fact, the assumption p to t (0) = p ® pE is quite unrealistic if the coupling is 
strong, as discussed below. Nonetheless, it certainly represents a valid initial 
state. 
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relevant at cu^t = 1, provided Xv&c (t) is then still appreciable, i.e. for 
4:\gk\ 2 /h 2 uj1 <C 1. Being a sum over periodic functions, Xvac(^) is quasi- 
periodic, that is, it will come back arbitrarily close to unity after a large 
period (which increases exponentially with the number of modes). These 
somewhat artificial Poincare recurrences vanish if we replace the sum over 
the discrete modes by an integral over a continuum with mode density [i, 

poo 

y;/H)— / du^M/M, (5.4i) 

k J o 

for any function /. This way the coupling constants gk get replaced by the 
spectral density of the environment, 

J(u) = 4|i(w) \g(uj)\ 2 . (5.42) 

This function characterizes the environment by telling how effective the cou- 
pling is at a certain frequency. 

Thermal State 

If the environment is in a thermal state with temperature T, 

,-H E /k B T 

tr (e~ HE / fc B T 



(Is) 

=/>ih 



the suppression factor read^l 



(5.43) 



k k \ 

(5.44) 

This factor can be separated into its vacuum component (|5.40|) and a ther- 
mal component, = e~ ^(''e - ■ F ">(*) ; with the following definitions of the 
vacuum and the thermal decay functions: 

iW*):=£4| gfc l 2 l -gt fct) , (5-45) 



5 This can be found in a small exercise by using the Baker-Hausdorff relation 
with exp (afc>* — a*b) = cxp(— |a| 2 /2) exp (ab') exp (— a*b), and the fact that 
coherent states satisfy the eigenvalue equation b\/3) = /3\0), have the number 
representation (n\/3) = exp(— \/3\ 2 /2)(3™ /y/n\, and form an over-complete set with 

I=TT- 1 fd 2 (}\P){P\- 
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5.2.2 The Continuum Limit 

Assuming that the field modes are sufficiently dense we replace their sum by 
an integration. Noting (|5.41|) . (15.421) we have 

F vac (t) — r dujJiu) IzggM (5.47) 
Jo n, U) 

t, , s f°° , r, s 1 - COs(wt) 

F th (t) — » j d^JH fi2 j ; 

So far, the treatment was exact. To continue we have to specify the spectral 
density in the continuum limit. A typical model takes g oc y/u, so that the 
spectral density of a d-dimensional field can be written as [17] 

J(w)=aw(— J e~ u ^ (5.49) 

with "damping strength" a > 0. Here, ui c is a characteristic frequency "cutoff" 
where the coupling decreases rapidly, such as the Debye frequency in the case 
of phonons. 



Ohmic Coupling 

For d = 1 the spectral density (|5.49| increases linearly at small w ("Ohmic 
coupling" ) . One finds 

Fvac(*) = ^l0g(l+C^ 2 ), (5.50) 

which bears a strong oj c dependence. Evaluating the second integral requires 
to assume that the cutoff lo c is large compared to the thermal energy, kT <C 
hhj c : 

Fa(t) ^ loe (^My „ 51) 

Here bp — h/ (nksT) is a thermal quantum time scale. The corresponding 
frequency u>i = 2/tx is called the (first) Matsubara frequency, which also 
shows up if imaginary time path integral techniques are used to treat the 
influence of bosonic field couplings [T7]. For large times the decay function 
F t h(i) shows the asymptotic behavior 

*th(*)~4^- [asi->oo]. (5.52) 

It follows that the decay of coherence is characterized by rather different 
regimes. In the short time regime (t < w" 1 ) we have the perturbative behavior 



(5.53) 
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which can also be obtained from the short-time expansion of the time evolu- 
tion operator. Note that the decay is here determined by the overall width tu c 
of the spectral density. The intermediate region, uj^ 1 < t < uj^ 1 , is dominated 
by F vac (t) and called the vacuum regime, 

F(t) ~ log (« c t) [for uj- 1 < t < ui 1 ] . 

Beyond that, for large times the decay is dominated by the thermal suppres- 
sion factor, 

[for^r 1 ^]. (5-54) 

In this thermal regime the decay shows the exponential behavior typical for 
the Markovian master equations discussed below. Note that the decay rate 
for this long time behavior is determined by the low frequency behavior of 
the spectral density, characterized by the damping strength a in (|5.49[) . and 
is proportional to the temperature T. 

Super-Ohmic Coupling 

For d = 3, the case of a "super-Ohmic" bath, the integrals (|5.4T[> . (|5.48|) can 
be calculated without approximation. We note only the long-time behavior 
of the decay, 

hm F(t) =2a( $ ( 1 + ^) < oo . (5.55) 
t^oo y fujjf. J \ nuj c J 

Here ip ( z ) stands for the Digamma function, the logarithmic derivative of 
the gamma function. Somewhat surprisingly, the coherences do not get com- 
pletely reduced as t — > oo, even at a finite temperature. This is due to the sup- 
pressed influence of the important low-frequency contributions to the spectral 
density in three dimensions (as compared to lower dimensions). While such 
a suppression of decoherence is plausible for intermediate times, the limiting 
behavior (|5.55|) is clearly a result of our simplified model assumptions. It will 
be absent if there is a small anharmonic coupling between the bath modes 
|18j or if there is a small admixture of different couplings to H; n t- 

Decoherence by "Vacuum Fluctuations"? 

The foregoing discussion seems to indicate that the "vacuum fluctuations" 
attributed to the quantized field modes are responsible for a general decoher- 
ence process, which occurs at short time scales even if the field is in its ground 
state. This ground state is non-degenerate and the only way to change it is 
to increase the energy of the field. But in our model the interaction Hamilto- 
nian if; nt commutes with the system Hamiltonian, so that it cannot describe 
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energy exchange between qubit and field. One would therefore expect that 
after the interaction the field has the same energy as before, so that an initial 
vacuum state remains unchanged and decoherence cannot take place. 

This puzzle is resolved by noting that the initial state pE — &) k |0)(0|fc 
is an eigenstate only in the absence of the coupling H[ nt , but not of the 
total Hamiltonian. By starting with the product state ptot(O) = p® pe we 
do not account for this possibly strong coupling. At an infinitesimally small 
time later, system and field thus suddenly feel that they are coupled to each 
other, which leads to a renormalization of their energies (as described by the 
Lamb shift discussed in Sect. I5.4.T]) . The factor Xvac in Eq. (|5.40|) describes 
the 'initial jolt' produced by this sudden switching on of the coupling. 

It follows that the above treatment of the short time dynamics, though 
formally correct, does not give a physically reasonable picture if the system 
state is prepared in the presence of the coupling. In this case, one should 
rather work with the eigenstates of the total Hamiltonian, often denoted as 
"dressed states" . If we start with a superposition of those two dressed states, 
which correspond in the limit of vanishing coupling to the two system states 
and the vacuum field, the resulting dynamics will show no further loss of 
coherence. This is consistent with the above notion that at zero temperature 
elastic processes cannot lead to decoherence [T5] . 

5.2.3 Dephasing of N Qubits 

Let us now discuss the generalization to the case of N qubits which do not 
interact directly among each other. Each qubit may have a different coupling 
to the bath modes. The system Hamiltonian is then 

N-i % N-x 

Htot = E + E + E E (^M + bPn*) . 

3=0 k j=Q k 

Ho 

(5.56) 

Similar to above, the time evolution in the interaction picture reads 

( N-l 
^E^E («J°bt _ [«y>]*b fc ) 
j=0 k 

where the displacement of the field modes now depends on the TV-qubit state. 

As an example, we take N = 2 qubits and only a single vacuum mode. 
For the inital qubit states 

|0) = c n | TTRcoolU} (5.57) 

and 




l^)=cio| TD + coil IT) 



(5.58) 
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we obtain, respectively, 

(t) 



U 

and 



= ch|TT)|- 



,( 2 ) 



cool 44)1- 



W (t) - a< 2 ) 



U|V)|0)e = c 1o |T4)|- 



,(D 



(*) 



,(2) 



coil 44)1- 



(*) 



,(2) 



where the a*- 1 ) (t) and a' 2 ) (t) are the field displacements (|5.34p due to the 
first and the second qubit. 

If the couplings to the environment are equal for both qubits, say, because 
they are all sitting in the same place and seeing the same field, we have 
aS 1 ' (t) — (t) = a (t). In this case, states of the form \<f>) are decohered 
once the factor (a(t) \ — a(t))v = exp(— 2 \a (t)\ 2 ) is approximately zero. 
States of the form on the other hand, are not affected at all, and one 
says that the {\ifi}} span a (two-dimensional) decoherence-free subspace. It 
shows up because the environment cannot tell the difference between the 
states | 44) an( i I 44) if h couples only to the sum of the excitations. 

For an arbitrary number of qubits, using an TV-digit binary notation, e.g. 
| T44) = |101 2 ) = |5), one has 



(m\~p{t)\n) 




■E(«£°(*) b I 



Pth 



(5.59) 



where mj € {0, 1} indicates the j-th digit in the binary representation of the 
number m. 

We can distinguish different limiting cases: 



Qubits Feel the Same Reservoir 

If the separation of the qubits is small compared to the wave lengths of the 

(i) 

field modes they are effectively interacting with the same reservoir, aY = a/.. 
One can push the j-summation to the a's in this case, so that, compared to 
the single qubit, one merely has to replace by XX m j — n j) a k ■ We find 



Xmn(t) = exp 



N-1 
3=0 



{Fy^{t) + Fth(*)) 



(5.60) 



with F vac (t) and F th (t) given by (|CT7|) and ![5"3g]) . 
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Hence, in the worst case, one observes an increase of the decay rate by 
iV 2 compared to the single qubit rate. This is the case for the coherence 
between the states |0) and 12^ — 1), which have the maximum difference 
in the number of excitations. On the other hand, the states with an equal 
number of excitations form a decoherence- free subspace in the present model, 
with a maximal dimension of (^ 2 ) ■ 

Qubits See Different Reservoirs 

In the other extreme, the qubits are so far apart from each other that each 
field mode couples only to a single qubit. This suggests a re- numbering of the 
field modes, 

a k ] — ► a k 3 , 

and leads, after transforming the j-summation into a tensor-product, to 

JV-l 

\mn(t) = J] tr ( ® (K ~ nj) a kj (*)) 

2 (^vacW + *th(t))) 

(FvacW+^thW)^). (5.61) 
Hamming distance 

Hence, the decay of coherence is the same for all pairs of states with the 
same Hamming distance. In the worst case, we have an increase by a factor 
of N compared to the single qubit case, and there are no decoherence-free 
subspaces. 

An intermediate case is obtained if the coupling depends on the position r j 
of the qubits. A reasonable model, corresponding to point scatterings of fields 
with wave vector fc, is given by gjp = g k exp (ifc • rj), and its implications 
are studied in [20] , 

The model for decoherence discussed in this section is rather exceptional 
in that the dynamics of the system can be calculated exactly for some choices 
of the environmental spectral density. In general, one has to resort to approx- 
imate descriptions for the dynamical effect of the environment; we turn to 
this problem in the following section. 




5.3 Markovian Dynamics of Open Quantum Systems 



Isolated systems evolve, in the Schrodinger picture and for the general case 
of mixed states, according to the von Neumann equation, 
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d t p=-[H,p}. (5.62) 

One would like to have a similar differential equation for the reduced dynam- 
ics of an "open" quantum system, which is in contact with its environment. 
If we extend the description to include the entire environment He and its 
coupling to the system, then the total state in Htot — H <S> He evolves uni- 
tarily. The partial trace over He gives the evolved system state, and its time 
derivative reads 

%>= — tr E (u tot (<) / 5tot(0)U t t ot W) = T^tr E ([H tot ,p tot ]) . (5.63) 

This exact equation is not closed and therefore not particularly helpful as it 
stands. However, it can be used as the starting point to derive approximate 
time evolution equations for p, in particular, if it is permissible to take the 
initial system state to be uncorrelated with the environment. 

These equations are often non-local in time, though, in agreement with 
causality, the change of the state at each point in time depends only on the 
state evolution in the past. In this case, the evolution equation is called a 
generalized master equation. It can be specified in terms of superoperator- 
functionals, i.e. linear operators which take the density operator p with its 
past time-evolution until time t and map it to the differential change of the 
operator at that time, 

d t p = K [{ Pt : T < t}} . (5.64) 

An interpretation of this dependence on the system's past is that the en- 
vironment has a memory, since it affects the system in a way which de- 
pends on the history of the system-environment interaction. One may hope 
that on a coarse-grained time-scale, which is large compared to the inter- 
environmental correlation times, these memory effects might become irrel- 
evant. In this case, a proper master equation might be appropriate, where 
the infinitesimal change of p depends only on the instantaneous system state, 
through a Liouville super-operator £, 

d t p = Cp . (5.65) 

Master equations of this type are also called Markovian, because of their 
resemblance to the differential Chapman-Kolmogorov equation for a classical 
Markov process. However, since a number of approximations are involved 
in their derivation, it is not clear whether the corresponding time evolution 
respects important properties of a quantum state, such as its positivity. We 
will see that these constraints restrict the possible form of C rather strongly. 

5.3.1 Quantum Dynamical Semigroups 

The notion a of quantum dynamical semigroup permits a rigorous formulation 
of the Markov assumption in quantum theory. To introduce it we first need a 
number of concepts from the theory of open quantum systems [TTJ [STJ H3] • 
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Dynamical Maps 

A dynamical map is a one-parameter family of trace-preserving, convex linear, 
and completely positive maps (CPM) 

W t :po^Pt, for<eR+, (5.66) 

satisfying Wo = id. As such, it yields the most general description of a time 
evolution which maps an arbitrary initial state po to valid states at later 
times. 

Specifically, the condition of trace preservation guarantees the normaliza- 
tion of the state, 

ti(p t ) = 1 , 

and the convex linearity, i.e. 

m (Xpo + (1 - A) p' ) = XWt (po) + (1 - A) W t (p' ) for all ^ A < 1 , 

ensures that the transformation of mixed states is consistent with the classical 
notion of ignorance. The final requirement of complete positivity is stronger 
than mere positivity of Wt (po)- It means that in addition all the tensor 
product extensions of Wt to spaces of higher dimension, defined with the 
identity map id ox t, are positive, 

W t <8> id cxt > , 

that is, the image of any positive operator in the higher dimensional space 
is again a positive operator. This guarantees that the system state remains 
positive even if it is the reduced part of a non-separable state evolving in a 
higher dimensional space. 

Kraus Representation 

Any dynamical map admits an operator-sum representation of the form (|5.8p 

N 

W t (p) = X>*(*)P W fc(*) (5.67) 
fe=i 

with the completeness relatior@ 

N 

^W£(i)W fc (i) =1. (5.68) 

k=l 

6 In case of a trace- decreasing, convex linear, completely positive map the condition 
(^gg)> is replaced by J2 k W it(*) W k(*) < I > i- e - thc operator I - E fc W I (*)W fc (t) 
must be positive. 
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The number of the required Kraus operators W/j(i) is limited by the dimen- 
sion of the system Hilbert space, N ^ dim (H) 2 (and confined to a count- 
able set in case of an infinite-dimensional, separable Hilbert space), but their 
choice is not unique. 

Semigroup Assumption 

We can now formulate the assumption that the {Wt '■ t € Rj} form a con- 
tinuous dynamical semigroup [211 123] : 

Wt 3 (W tl (-))=Wt l+ t a (-) forallti,t2>0 (5.69) 

and Wo = id. This statement is rather strong, and it is certainly violated 
for truly microscopic times. But it seems not unreasonable on the level of 
a coarse-grained time scale, which is long compared to the time it takes for 
the environment to "forget" the past interactions with the system due to the 
dispersion of correlations into the many environmental degrees of freedom. 

For a given dynamical semigroup there exists, under rather weak condi- 
tions, a generator, i.e. a superoperator £ satisfying 

W t = e ct for t > . (5.70) 

In this case Wt(p) is the formal solution of the Markovian master equation 

ssm . 

Dual Maps 

So far we used the Schrodingcr picture, i.e. the notion that the state of an open 
quantum system evolves in time, p t = Wt(po)- Like in the description of closed 
quantum systems, one can also take the Heisenberg point of view, where the 
state does not evolve, while the operators A describing observables acquire 
a time dependence. The corresponding map VVj : Ao i— ► At is called the dual 
map, and it is related to W* by the requirement tr(AWt(p)) = tr(/oWf(A)). 
In case of a dynamical semigroup, Wf = exp (£"t) , the equation of motion 
takes the form 9jA = £"A, with the dual Liouville operator determined by 
tr(A£ (p)) = tr(/o£" (A)). From a mathematical point of view, the Heisenberg 
picture is much more convenient since the observables form an algebra, and 
it is therefore preferred in the mathematical literature. 

5.3.2 The Lindblad Form 

We can now derive the general form of the generator of a dynamical semi- 
group, taking dim (Ti) = d < oo for simplicity The bounded operators 

7 The inverse element required for a group structure is missing for general, irre- 
versible CPMs. 
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on 7i then form a cP-dimensional vector space which turns into a Hilbert 
space, if equipped with the Hilbert-Schmidt scalar product (A, B) := tr(A^ B). 
Given an orthonormal basis of operators {E^ : 1 ^ j ^ d 2 } C L (H), 

(Ei,E i ):=tr(EjE J ) = ^, (5.71) 

any Hilbert-Schmidt operator Wfe can be expanded as 

d 2 

W fe = J](E„W fe )E J ;. (5.72) 

3=1 

We can choose one of the basis operators, say the d 2 -th, to be proportional 
to the identity operator, 

E d2 = -J=I , (5.73) 

so that all other basis elements are traceless, 

fo for j = 1,...,(P - 1 , 

lt(E ' ) = {^ Jj-f. (5 ' 74) 

Representing the superoperator of the dynamical map (|5.67fl in the {Ej} basis 
we have 

d 2 

W t (p) = c «(*)EipE] (5.75) 

i ,3 = 1 

with a time dependent, hermitian and positive coefficient matrix, 

N 

*i(t) = E( E ^ w *(*))(Ej, W fc (i))* (5.76) 

k=l 

(positivity can be checked in a small calculation). We can now calculate the 
semigroup generator in terms of the differential quotient by writing the terms 
including the element E d 2 separately: 
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Lp = lim 

r^O T 



\c d 2 d 2{T) - 1 t ,.„ v^ c j«p( t )c - , V^C,Pj( t )i 



= lim 



-p+hm > — Ejp + phm > — j= — E\ 



Bei(W) 



BteL(«) 



V lim^E^E] 



i ,3=1 " 



= c p + Bp + pBt + a H E ipEj 



i.j=i 



^ [H , p] + 1 (Gp + pG) + ^ aij E iP Et . 



(5.77) 



In the last equality the following hermitian operators with the dimension of 
an energy were introduced: 

G = ^(B + Bt+ Co ), 

H = £(B-B+) . 
2i v ; 

By observing that the conservation of the trace implies tr(£p) = 0, one can 
relate the operator G to the matrix a = (aij), since 



= tr(£p) = + tr 
must hold for all p. It follows that 



2G V? ptp 



This leads to the first standard form for the generator of a dynamical semi- 
group: 

d 2 -l 

C P = ^[H,p] + E a « ( E *<° E ] - Ie+E^p- IpE+E*) . (5.78) 



unitary part 



incoherent part 



The complex coefficients have dimensions of frequency and constitute a 
positive matrix ex. 
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The second standard form or Lindblad form is obtained by diagonaliz- 
ing the coefficient matrix a. The corresponding unitary matrix U satisfying 
UaU^ = diag(7i, . . . ,7d2_i) allows to define the dimensionless operators 

Lfe := Yfj=i Ej^jfc 80 tnat E .j = Yfk=i \-kU kj and therefor^] 



£p--[H,p] + ^ 7fe 
fc=i 




(5.79) 



with N d 2 — 1. This shows that the general form of a generator of a 
dynamical semigroup is specified by a single hermitian operator H, which is 
not necessarily equal to the Hamiltonian of the isolated system, see below, 
and at most d 2 — 1 arbitrary operators with attributed positive rates 7&. 
These are called Lindblad operator^ or jump operators, a name motivated 
in the following section. 

It is important to note that a given generator C does not determine the 
jump operators uniquely. In fact, the equation is invariant under the trans- 
formation 

U -» L fc + c k , (5.80) 

h ^ h +|e^- ( c i L i- c i L i) . ( 5 - si ) 

3 

with Ck £ C, so that the can be chosen traceless. In this case, the only 
remaining freedom is a unitary mixing, 

^HU^Y, U 'i3^3- (5-82) 

3 

If C shows an additional invariance, e.g. with respect to rotations or trans- 
lations, the form of the Lindblad operators will be further restricted, see 
e.g. |23| . 



5.3.3 Quantum Trajectories 

Generally, if we write the Liouville superoperator C as the sum of two parts, 
Co and iS, then the formal solution (|5.70| of the master equation (|5.65|) can 
be expressed as 

8 It is easy to see that the dual Liouville operator discussed in Sect. 15.3.11 reads, 
in Lindblad form, 

1 - / 1 1 \ 

r 8 (A) = -[A, H] + £ 7fc 4AL* - -4UA - -ALlU . 

Note that this implies £ B (I) = 0, while C (I) = £ fc 7fc [U, L.£], and tr (£X) = 0. 

9 Lindblad showed in 1976 that the form (15.79[) is obtained even for infinite- 
dimensional systems provided the generator C is bounded (which is usually not 
the case). 
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E 
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ra=l 



2 / dt„ / dt„_! • • • / d«i 



/0 JO Jo 

xe A^t-tn)£ e A)(t»-t»-i)s... e £o(*2-ti) < s e £oti _ ( 5 83 ) 



IS 



The step from the second to the third line, where t n+ ^ kj j(n + Ej fcj) 
replaced by n time integrals, can be checked by induction. 

The form (|5.83|) is a generalized Dyson expansion, and the comparison 
with the Dyson series for unitary evolutions suggests to view exp (Cot) as 
an "unperturbed" evolution and S as a "perturbation" , such that the exact 
time evolution Wt is obtained by an integration over all iterations of the 
perturbation, separated by the unperturbed evolutions. 

The particular Lindblad form (|5.79[) of the generator suggests to introduce 
the completely positive jump superoperators 

C kP = TfeUpl-I , (5.84) 
along with the non-hermitian operator 

•fi N 

Hc = H- yX>4L fc . (5.85) 

fc=i 

The latter has a negative imaginary part, Im (H<c) < 0, and can be used to 
construct 

C p=^(H c p-pHl) . (5.86) 
It follows that the sum of these superoperators yields the Liouville operator 

(529 

N 

£= £ + J2C k . (5.87) 

fc=i 
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Of course, neither £ nor <5 = SfcLi generate a dynamical semigroup. 
Nonetheless, they are useful since the interpretation of (|5.83f> can now be 
taken one step further. We can take the point of view that the Ck with k ^ 1 
describe elementary transformation events due to the environment ( "jumps" ) , 
which occur at random times with a rate 7&. A particular realization of n 
such events is specified by a sequence of the form 

Rn = {ti, ki;t 2 , fc 2 ; . . . .; t n , k n ) . (5.88) 

The attributed times satisfy < ti < < t n < t, and the kj € {1, . . . , N} 

indicate which kind of event "took place" . We call a record of length n. 

The general time evolution W t can thus be written as an integration over 
all possible realizations of the jumps, with the "free" evolution exp (Cqt) in 
between, 



n=l ' 
X 

{Rn} 



Wt = e Cat + V / dt n / dt„_x • ■ ■ / d^ 
=1 Jo Jo Jo 

J2 e £ °( t -*")£ fc „e £ °( t ''-*"- 1 )A n _ 1 • • -e^-^Ck^ 



(5.89) 

As a result of the negative imaginary part in (|5.85p the exp (Co T ) are trace 
decreasincP^ completely positive maps, 

e £or p = exp^-^H c ^pexp^H^ >0, (5.90) 
d N 

— tr (e CoT p) = tr (C e c ° T p) = - ^ tr(£fce^p) < . (5.91) 



fc-i >0 



It is now natural to interpret tr (e £o *p) as the probability that no jump occurs 
during the time interval t, 

Prob(i?^|p) := tr(e £ot p) • (5.92) 

To see that this makes sense, we attribute to each record of length n a rt-time 
probability density. For a given record R„ we define 

prob«|p) :=tr(/Cfl*p) , (5.93) 

in terms of the superoperators from the second line in (|5.89[) , 

JC Rt :=e £ °( t - t ")/: fc , i e £t, (*"-*"- 1 )/: fc „_ 1 ---e £o ( t2 - tl )/: fcl e £otl . (5.94) 



See the note [6] 
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This is reasonable since the ICr^ are completely positive maps that do not 
preserve the trace. Indeed, the probability density for a record is thus de- 
termined both by the corresponding jump operators, which involve the rates 
7fc, and by the e c ° T p, which account for the fact that the likelihood for the 
absence of a jump decreases with the length of the time interval. 

This notion of probabilities is consistent, as can be seen by adding the 
probability (|5.92|) for no jump to occur during the interval (0; t) to the inte- 
gral over the probability densities (|5.93|) of all possible jump sequences. As 
required, the result is unity, 



00 r* pin rt°2 

Prob(i^|p)+y" / dt n / dt n _i--- / dii^prob« 

JO JO JO y Rt y 



for all p and t ^ 0. This follows immediately from the trace preservation of 
the map 1(5155)) . 

It is now natural to normalize the transformation defined by the ICr^ ■ 
Formally, this yields the state transformation conditioned to a certain record 
i?^. It is called a quantum trajectory, 

T^):,-^. (5.95) 

Note that this definition comprises the trajectory corresponding to a null- 
record R , where IC R t = exp(£oi)- These completely positive, trace-pre- 
serving, nonlinear maps p 1— > T (p\R n ) are defined for all states p that yield 
a finite probability (density) for the given record i?^, i.e. if the denominator 
in (|5.95p does not vanish. 

Using these notions, the exact solution of a general Lindblad master equa- 
tion (|5.83p may thus be rewritten in the form 

p t = Prob(^) T(p\R ) 

00 pt pt n pt-2 

+ dt n dtn-i-- dhJ2 prob (R n \p) T(p\R n ) . 

n=l J JO JO {Rn} 

(5.96) 

It shows that the general Markovian quantum dynamics can be understood 
as a summation over all quantum trajectories T (fl\R n ) weighted by their 
probability (density). This is called a stochastic unravelling of the master 
equation. The set of trajectories and their weights are labeled by the possi- 
ble records (|5.88[) . and determined by the Lindblad operators of the master 
equation (|5.79[) . 

The semigroup property described by the master equation shows up if a 
record is formed by joining the records of adjoining time intervals, (0; t') and 

(f;t), 
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- (^° ;t 'U ( n ■ (5-97) 
As one expects, the probabilities and trajectories satisfy 

W ob(R^l\p) = prob(i?I° ;t,) |p) prob(i?^' ;t) |T(p|^i 0; *' ) )) , (5-98) 

and 

T(p\R ( n+l) = TiTiplR^'^U^h . (5.99) 

Note finally, that the concept of quantum trajectories fits seamlessly into 
the framework of generalized measurements discussed in Sect. [BTTT51 In par- 
ticular, the conditioned state transformation T (p\R^) has the form (|5 . 18[) of 
an efficient measurement transformation, 

M R t oML 

T (p\Ri ) = — 2»_ (5.100) 

tr \M r Rl M K p) 

with compound measurement operators 

Mk .- e -iH«c(t-t„)/fi Lfcra . . . Lfc2e -iH c (t 2 - tl )/ft Lfcie -iH ctl /ft _ ( 5AQ l) 

This shows that we can legitimately view the open quantum dynamics gen- 
erated by C as due to the continuous monitoring of the system by the en- 
vironment. We just have to identify the (aptly named) record i?^ with the 
total outcome of a hypothetical, continuous measurement during the interval 
(0;t). The jump operators Lk then describe the effects of the corresponding 
elementary measurement event£3 ("clicks of counter fc"). Since the absence 
of any click during the "waiting time" r may also confer information about 
the system, this lack of an event constitutes a measurement as well, which is 
described by the non-unitary operators exp (— iH<cr/fi). A hypothetical de- 
mon, who has the full record R^ available, would then be able to predict the 
final state T (p\R^). In the absence of this information we have to resort to 
the probabilistic description (|5.96[) weighting each quantum trajectory with 
its (Bayesian) probability. 

We can thus conclude that the dynamics of open quantum dynamics, and 
therefore decoherence, can in principle be understood in terms of an infor- 
mation transfer to the environment. Apart from this conceptual insight, the 
unravelling of a master equation provides also an efficient stochastic simula- 
tion method for its numerical integration. In these quantum jump approaches 

11 For all these appealing notions, it should be kept in mind that the Ck are not 
uniquely specified by a given generator C, see (I5.80|) - (|5.82|l . Different choices 
of the Lindblad operators lead to different unravellings of the master equation, 
so that these hypothetical measurement events must not be viewed as "real" 
processes. 
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[551 123 , which are based on the observation that the quantum trajectory 
(1 5 . 9 5 [) of a pure state remains pure, one generates a finite ensemble of trajec- 
tories such that the ensemble mean approximates the solution of the master 
equation. 

5.3.4 Exemplary Master Equations 

Let us take a look at a number of very simple Markovian master equations^, 
which are characterized by a single Lindblad operator L (together with a her- 
mitian operator H). The first example gives a general description of dephas- 
ing, while the others are empirically known to describe dissipative phenomena 
realistically. We may then ask what they predict about decoherence. 

Dephasing 

The simplest choice is to take the Lindblad operator to be proportional to 
the Hamiltonian of a discrete quantum system, i.e. to the generator of the 
unitary dynamics, L = ^/t~H. The Lindblad equation 



is immediately solved in the energy eigenbasis, H = J2m E m \rn){m\ : 



As we expect from the discussion of qubit dephasing in Sect. [521 the energy 
eigenstates are unaffected by the non-unitary dynamics if the environmental 
effect commutes with H. The coherences show the exponential decay that we 
found in the "thermal regime" (of times t which are long compared to the 
inverse Matsubara frequency). The comparison with (|5.54p indicates that 7 
should be proportional to the temperature of the environment. 

Amplitude Damping of the Harmonic Oscillator 

Next, we choose H to be the Hamiltonian of a harmonic oscillator, H = fkua^a, 
and take as Lindblad operator the ladder operator, L = a. The resulting 
Lindblad equation is known empirically to describe the quantum dynamics 
of a damped harmonic oscillator. 



See also Sect. 6 in Cord Muller's contribution for a discussion of the master 
equation describing spin relaxation. 

As an exception, 7 does not have the dimensions of a rate here (to avoid clumsy 
notation). 




(5.102) 




(5.103) 
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Choosing as initial state a coherent state, see (|5.37p and (|5 . 1 14)) below, 

Po = \a ){a \ = D(ao)|0)<0|Dt(ao) 

= e- |ao|2 exp(a o a t )|0)(0|exp(a5a) , (5.104) 

we are faced with the exceptional fact that the state remains pure during the 
Lindblad time evolution. Indeed, the solution of the Lindblad equation reads 

pt = \a t )(a t \ (5.105) 

with 



at 



ag exp (^—iuit — ~ij , (5.106) 



as can be verified easily using (|5. 104|) . It describes how the coherent state 
spirals in phase space towards the origin, approaching the ground state as 
t — > co. The rate 7 is the dissipation rate since it quantifies the energy loss, 
as shown by the time dependence of the energy expectation value, 

(a t |H|a*) = e-T rt (ao|H|ao>. (5.107) 
Superposition of Coherent States 

What happens if we start out with a superposition of coherent states, 

\ 1 fo) = -L-(\a a ) + \0 o )) (5.108) 



with J\f = 2 + 2Re(ao|/3o), in particular, if the separation in phase space is 
large compared to the quantum uncertainties, | cto — /3o | 3> 1 ? The initial 
density operator corresponding to (|5. 108(1 reads 

po = ^ (\a )(a \ + \p )(po\ + co\a )(/3 \ + c*\p )(a \) (5.109) 

with Co = 1. One finds that the ansatz 

Pt = jf(\at)(at\ + \fr)(fr\+c t \c h )(fr\ (5.110) 

solves the Lindblad equation with (|5.106[) . provided 



c t = c exp 



"2 K - A) I 2 + ilm(a $)) 



(1-e-T*)] . (5.111) 



That is, while the coherent "basis" states have the same time evolution as 
in (|5.105[) . the initial coherence Co gets additionally suppressed. For times 
that are short compared to the dissipative time scale, t 7 _1 , we have an 
exponential decay 
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lotl = |co|cxp(- ||a -/3 |V) (5.112) 

7doco 

with a rate 7deco- For macroscopically distinct superpositions, where the 
phase space distance of the quantum states is much larger than their un- 
certainties, \cto — A) I ^ 1: the decoherence rate 7doco can be much greater 
than the dissipation rate, 

— ^K-/?o| 2 »l. (5.113) 
7 2 

This quadratic increase of the decoherence rate with the separation between 
the coherent states has been confirmed experimentally in a series of beautiful 
cavity QED experiments in Paris, using field states with an average of 5-9 
photons [251125]. 

Given this empiric support we can ask about the prediction for a material, 
macroscopic oscillator. As an example, we take a pendulum with a mass of 
m =100 g and a period of 2tt/uj = Is, and assume that we can prepare it in 
a superposition of coherent states with a separation of x =lcm. The mode 
variable a is related to position and momentum by 




(5.114) 



so that we get the prediction 

7dcco ~ 10 30 7 . 

This purports that even with an oscillator of enormously low friction cor- 
responding to a dissipation rate of 7 =l/year the coherence is lost on a 
timescale of 10 -22 s - in which light travels the distance of about a nuclear 
diameter. 

This observation is often evoked to explain the absence of macroscopic su- 
perpositions. However, it seems unreasonable to assume that anything phys- 
ically relevant takes place on a timescale at which a signal travels at most by 
the diameter of an atomic nucleus. Rather, one expects that the decoherence 
rate should saturate at a finite value if one increases the phase space distance 
between the superposed states. 



Quantum Brownian Motion 

Next, let us consider a particle in one dimension. A possible choice for the 
Lindblad operator is a linear combination of its position and momentum 
operators, 
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Here p t h is a momentum scale, which will be related to the temperature of the 
environment below. The hermitian operator is taken to be the Hamiltonian of 
a particle in a potential V (x), plus a term due to the environmental coupling, 

2 

H = ^ + V(x) + |(xp + px) . (5.116) 

This additional term will be justified by the fact that the resulting Lindblad 
equation is almost equal to the Caldeira-Leggett master equation. The latter 
is the high-temperature limit of the exact evolution equation following from 
a harmonic bath model of the environment [301 US], see Sect. 15.4.11 It is 
empirically known to describe the frictional quantum dynamics of a Brownian 
particle, and, in particular, for t — > oo it leads to the canonical Gibbs state 
in case of quadratic potentials. 

The choices (15 . 1 15[) and (|5.116[) yield the following Lindblad equation: 

Caldcira— Lcggctt master equation 



dtPt = h%- + V(x), Pt ] + 1 [x, PPt + p tP ]-?if [x, [x, Pt] 
in, 2m m 2 nr 



dissipation position localization 

-I-L[P,[P,Pt]] ■ (5.H7) 

Z Pth 

The three terms in the upper line (with p t h from (|5.123| ) constitute the 
Caldeira-Leggett master equation. It is a Markovian, but not a completely 
positive master equation. In a sense, the last term in (15.117|) adds the minimal 
modification required to bring the Caldeira-Leggett master equation into 
Lindblad form [2 [32]. 

To see the most important properties of (|5.117p let us take a look at 
the time evolution of the relevant observables in the Heisenberg picture. As 
discussed in Sect. 15.3.11 the Heisenberg equations of motion are determined 
by the dual Liouville operator C.K In the present case, it takes the form 

(A) = i [AX + V(x)} - 1 (p [x, A] + [x, A] p) - 1% [x, [x, A]] 
in 2m in 2 nr 



ll 



2 L P,[P,A]] . (5.118) 



It is now easy to see that 

r« (x 



_p_ 

m ' 

C) (p) = -V'(x) - 2 7 p . (5.119) 



Hence, the force arising from the potential is complemented by a frictional 
force which will drive the particle into thermal equilibrium. The fact that 
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this frictional component stems from the second term in (|5.117p indicates 
that the latter describes the dissipative effect of the environment. 

In the absence of an external potential, V = 0, the time evolution deter- 
mined by (|5.119p is easily obtained, since (£") (p) = (—27)" p for n € N: 



n=0 



p t =e~"p= > 1 — p = e 'p 



x t = e £ * t x = x+-^- T (>C»)- *(p 



1 *" / \ , P - Pi 



[for V = 0] 



m ^— ' n! 27m 

n=l ' 



(5.120) 

Note that, unlike in closed systems, the Heisenberg operators do not retain 
their commutator, [x t , p t ] ^ ih for t > (since the map Wf = exp (£*t) is 
non- unitary). Similarly, (p 2 )t 7^ (Pi) 2 for i > 0, so that the kinetic energy 
operator T = p 2 /2m has to be calculated separately. Noting 

£ tt (T) = - 4 7 T [for V = 0] , (5.121) 

we find 

T t = Ek + (j -EL\ e -^t [forV = 0l. (5.122) 
8m \ 8m / 

This shows how the kinetic energy approaches a constant determined by the 
momentum scale pth- We can now relate pth to a temperature by equating the 
stationary expectation value tr (pToo) = p 2 h /8m with the the average kinetic 
energy ^k^T in a one-dimensional thermal distribution. This leads to 



Pth 



2^/mkgT (5.123) 



for the momentum scale in (|5. 1 15|) . Usually, one is not able to state the 
operator evolution in closed form. In those cases it may be helpful to take 
a look at the Ehrenfest equations for their expectation values. For example, 
given (p 2 )t = 2m(T) t , the other second moments, (x 2 ) t and (px + xp) t form a 
closed set of differential equations. Their solutions, given in [2], yield the time 
evolution of the position variance <r 2 (t) = (x 2 ) t — (x) 2 . It has the asymptotic 
form 

al (t) ~ -^-t as t -> 00 , (5.124) 
which shows the diffusive behavior expected of a (classical) Brownian parti- 

ckB 

Let us finally take a closer look at the physical meaning of the third term 
in (|5.117[) . which is dominant if the state is in a superposition of spatially 



Note that the definition of 7 differs by a factor of 2 in part of the literature. 
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separated states. Back in the Schrodinger picture we have in position repre- 
sentation, p t (x,x') = (x\p t \x'}, 

2 

d t pt(x,x') = --^-(x-x') 2 pt(x,x') + [the other terms] . (5.125) 
2 7r 



The "diagonal elements" p(x,x) are unaffected by this term, so that it leaves 
the particle density invariant. The coherences in position representation, on 
the other hand, get exponentially suppressed, 

p t {x 7 x 1 ) = cxp (-7dccoi) Po{x, x 1 ) . (5.126) 

Again the decoherence rate is determined by the square of the relevant dis- 
tance \x — x'\, 

l^L = 4n {x ~/ )2 . (5.127) 
7 A\ 

Like in Sect. 15.3.41 the rate 7doco will be much larger than the dissipative 
rate provided the distance is large on the quantum scale, here given by the 
thermal de Broglie wavelength 

In particular, one finds 7d CC o 3^ 7 if the separation is truly macroscopic. 
Again, it seems unphysical that the decoherence rate does not saturate as 
\x — x'\ — > oo, but grows above all bounds. One might conclude from this 
that non-Markovian master equations are more appropriate on these short 
timescales. However, I will argue that (unless the environment has very spe- 
cial properties) Markovian master equations are well suited to study deco- 
herence processes, provided they involve an appropriate description of the 
microscopic dynamics. 



5.4 Microscopic Derivations 

In this section we discuss two important and rather different strategies to 
obtain Markovian master equations based on microscopic considerations. 

5.4.1 The Weak Coupling Formulation 

The most widely used form of incorporating the environment is the weak 
coupling approach. Here one assumes that the total Hamiltonian is "known" 
microscopically, usually in terms of a simplified model, 



Htot = H + He + Hi n t 
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and takes the interaction part H; nt to be "weak" so that a perturbative treat- 
ment of the interaction is permissible. 

The main assumption, called the Born approximation, states that H; nt 
is sufficiently small so that we can take the total state as factorized, both 
initially, ptot(O) = p(0) <£> pe, and also at t > in those terms which involve 
Hi n t to second order. 

Assumption 1 : ptot(i) — p(t) <8> Pe [to 2nd order in H; nt ] . (5.129) 

Here pe is the stationary state of the environment, [He,Pe] — 0. Like above, 
the use of the interaction picture is indicated with a tilde, cf. (|5 . 28[) . so that 
the von Neumann equation for the total system reads 

dtptot = TT[Hint(*),Ptat(*)] 

= ^[hnt{t),Ptot(0)] + -^J ds [Hint W, [Hint (s),/5 tot (s)]] . 

(5.130) 

In the second equation, which is still exact, the von Neumann equation in 
its integral version was inserted into the differential equation version. Using 
a basis of Hilbert-Schmidt operators of the product Hilbert space, see Sect. 
15.3.21 one can decompose the general Hi nt into the form 

HintW =^A fe (t)® B fe (i) (5.131) 

fc 

with Afc = A' k , Bfc = B^. The first approximation is now to replace /?tot(s) by 
p(s) ® Pe in the double commutator of (|5.130[) . where H; n t appears to second 
order. Performing the trace over the environment one gets 

d t p{t) = tr E (9t/5 to t) 

4E(B*(t)> m [A*(i),P'(0)] 
fc 

+ E{/ ds (B k (t)B e Xs)) PE {Mt)Ms)p(s) - A e ( s )p( s )A k (t)} 

C k i(t-s) 

+ h.c.|. (5.132) 

All the relevant properties of the environment are now expressed in terms of 
the (complex) bath correlation functions Ckt(t — s). Since [He,Pe] = 0, they 
depend only on the time difference t — s, 

C u {t) = tr (e iHET B fc e- iHET B^ E ) = (e iHET B fcG - iH ^B,) . (5.133) 
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This function is determined by the environmental state alone, and it is typi- 
cally appreciable only for a small range of r around r = 0. 

Equation (|5. 132|) has the closed form of a generalized master equation, 
but it is non-local in time, i.e. non-Markovian. Viewing the second term as a 
superoperator JC, which depends essentially ont-s we have 

W) = i[(H int (% E) p(0)]+ / dsK(t - s)p(s) , (5.134) 

disregarded 

where JC is a superoperator memory kernel of the form (|5.64|) . We may disre- 
gard the first term since the model Hamiltonian He can always be reformu- 
lated such that (Bk(t)) PE — 0. 

A naive application of second order of perturbation theory would now 
replace p(s) by the initial p(0). However, since the memory kernel is dominant 
at the origin it is much more reasonable to replace p(s) by p(t). The resulting 
master equation is local in time, 

d t p{t) S + dstC(t - s)\ p{t) . (5.135) 

It is called the Redfield equation and it is not Markovian, because the inte- 
grated superoperator still depends on time. Since the kernel is appreciable 
only at the origin it is reasonable to replace t in the upper integration limit 
by oo. 

These steps are summarized by the Born-Markov approximation: 

C-t i* CO 

Assumption 2 : / dsJC(t — s)p (s) = / dsJC(s)p(t) . (5.136) 
Jo Jo 

It leads from (|5. 134(1 to a Markovian master equation provided (Hi nt ) pE = 0. 

However, by no means is such a master equation guaranteed to be com- 
pletely positive. An example is the Caldeira-Leggett master equation dis- 
cussed in Sect. 15.3.41 It can be derived by taking the environment to be a 
bath of bosonic field modes whose field amplitude is coupled linearly to the 
particle's position operator. A model assumption on the spectral density of 
the coupling then leads to the frictional behavior of (I5.119j) [30l fT7] . 

A completely positive master equation can be obtained by a further sim- 
plification, the "secular" approximation, which is applicable if the system 
Hamiltonian H has a discrete, non-degenerate spectrum. The system opera- 
tors Afc can then be decomposed in the system energy eigenbasis. Combining 
the contributions with equal energy differences 

A fc ( W )= ]T (E\A k \E')\E){E'\ =AlH, (5.137) 

E'-E=hu> 



we have 
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A fe = ^A fe ( w ). (5.138) 

The time dependence of the operators in the interaction picture is now par- 
ticularly simple, 

Mt) =5>- iwt A fc H . (5.139) 

OJ 

Inserting this decomposition we find 

(5.140) 

with 

r w (w) = ^ / d S e iws (B fe ( s )B,(0)) PE . (5.141) 



ft 2 7 

For times t which are large compared to the time scale given by the smallest 
system energy spacings it is reasonable to expect that only equal pairs of 
frequencies to, u) contribute appreciably to the sum in (|5.140p . since all other 
contributions are averaged out by the wildly oscillating phase factor. This 
constitutes the rotating wave approximation, our third assumption 

Assumption 3 : ^ ^"-"^ftu, u') ~ ^ f(u, u) . (5.142) 
It is now useful to rewrite 

= \lkl{u) + iSkt(u) (5.143) 

with 7fc^(w) given by the full Fourier transform of the bath correlation func- 
tion, 

1m {lo) = r M (w) + r/ fc (w) = ^ J dte iuJt (§ fc (t)B/(Q)^ , (5.144) 
and the hermitian matrix Ski(to) defined by 

S M (w) = |(r w (a;) - f^w)) . (5.145) 

The matrix 7/^(0;) is pos^toJ*"^! so that we end up with a master equation of 
the first Lindblad form (|5.78|) . 

9tp(t) = - [HL»»nb,p(i)] +Z)7«(w)( Mw)p(t)A£(oO 



^A+(a;)A<(w)p(t) - ip(t)A£(u;)A«(«) 



(5.147) 



To see that the matrix (7 (w)) fc £ = 7w(o;) is positive we write 
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The hermitian operator 



fi$^Sjw(w)A£(a;)A«(w) 



(5.148) 



describes a renormalization of the system energies due to the coupling with 
the environment, the Lamb shift. Indeed, one finds [H, HLamb] = 0. 

Reviewing the three approximations (|5.129|) . (|5.136|) , (|5.142j) in view of 
the decoherence problem one comes to the conclusion that they all seem to 
be well justified if the environment is generic and the coupling is sufficiently 
weak. Hence, the master equation should be alright for times beyond the 
short-time transient which is introduced due to the choice of a product state 
as initial state. Evidently, the problem of non-saturating decoherence rates 
encountered in Sect. 15.3.41 is rather due to the linear coupling assumption, 
corresponding to a "dipole approximation" , which is clearly invalid once the 
system states are separated by a larger distance than the wavelength of the 
environmental field modes. 

This shows the need to incorporate realistic, nonlinear environmental cou- 
plings with a finite range. A convenient way of deriving such master equations 
is discussed in the next section. 

5.4.2 The Monitoring Approach 

The following method to derive microscopic master equations differs consid- 
erably from the weak coupling treatment discussed above. It is not based 
on postulating an approximate "total" Hamiltonian of system plus environ- 
ment, but on two operators, which can be characterized individually in an 
operational sense. This permits to describe the environmental coupling in 
a non-perturbative fashion, and to incorporate the Markov assumption right 
from the beginning, rather than introducing it in the course of the calculation. 



with C := h^Y.t «/B/(0). One can now check that due to its particular form 
the correlation function 



appearing in (|5.146p is of positive type, meaning that the n x n- matrices 

(/ (ti ~ • defined by an arbitrary choice of ti, , t n and n £ N are positive. 

According to Bochner's theorem [35] the Fourier transform of a function which 
is of positive type is positive, which proves the positivity of (|5.146|) . 



0,7«) = ^2v* k j ke (uj)v e 




(5.146) 
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The approach may be motivated by the observation made in Sects. 15.1.31 
and 15.3.31 that environmental decoherence can be understood as due to the 
information transfer from the system to the environment occurring in a se- 
quence of indirect measurements. In accordance with this, we will picture 
the environment as monitoring the system continuously by sending probe 
particles which scatter off the system at random times. This notion will be 
applicable whenever the interaction with the environment can reasonably 
be described in terms of individual interaction events or "collisions" , and it 
suggests a formulation in terms of scattering theory, like in Sect. 15.1.21 The 
Markov assumption is then easily incorporated by disregarding the change of 
the environmental state after each collision [34] . 

When setting up a differential equation, one would like to write the tem- 
poral change of the system as the rate of collisions multiplied by the state 
transformation due to an individual scattering. However, in general not only 
the transformed state will depend on the original system state, but also the 
collision rate, so that such a naive ansatz would yield a nonlinear equation, 
violating the basic principles of quantum mechanics. To account for this state 
dependence of the collision rate in a proper way we will apply the concept 
of generalized measurements discussed in Sect. 15.1.31 Specifically, we shall 
assume that the system is surrounded by a hypothetical, minimally invasive 
detector, which tells at any instant whether a probe particle has passed by 
and is going to scatter off the system, see Fig. 15.11 




Fig. 5.1. (a) In the monitoring approach the system is taken to interact at most 
with one environmental (quasi-)particle at a time, so that three-body collisions 
are excluded. Moreover, in agreement with the Markov assumption, it is assumed 
that the environmental particles disperse their correlation with the system before 
scattering again, (b) In order to consistently incorporate the state-dependence of 
the collision rate into the dynamic description of the scattering process, we imagine 
that the system is monitored continuously by a transit detector, which tells at a 
temporal resolution At whether a particle is going to scatter off the system, or not. 
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The rate of collisions is then described by a positive operator T acting in 
the system- probe Hilbert space. Given the uncorrelated state g tot = p® p%, 
it determines the probability of a collision to occur in a small time interval 
At, 

Prob {C At \p ® p E ) = Attr(r[p®p E ]) . (5.149) 

Here, pe is the stationary reduced single particle state of the environment. 
The microscopic definition of T will in general involve the current density 
operator of the relative motion and a total scattering cross section, see below. 

The important point to note is that the information that a collision will 
take place changes our knowledge about the state, as described by the gen- 
eralized measurement transformation (|5. 16|1 . At the same time, we have to 
keep in mind that the measurement is not real, but is introduced here only 
for enabling us to account for the state dependence of the collision probabil- 
ity. It is therefore reasonable to take the detection process as efficient, see 
Sect. 15.1.31 and minimally-invasive, i.e. U a = I in (|5.20p . so that neither 
unnecessary uncertainty nor a reversible back-action is introduced. This im- 
plies that after a (hypothetical) detector click, but prior to scattering, the 
system-probe state will have the form 

ri/2 0f .rV2 

M (ftatlCW) = 7 ; r v • (5.150) 
tr (Iptot) 

This measurement transformation reflects our improved knowledge about the 
incoming two-particle wave packet, and it may be viewed as enhancing those 
parts which are heading towards a collision. Similarly, the absence of a detec- 
tion event during At changes the state, and this occurs with the probability 
Prob (C At ) = 1 - Prob (C At ). 

Using the state transformation (|5.150p we can now formulate the uncon- 
ditioned system-probe state after a coarse grained time At as the mixture 
of the colliding state transformed by the S-matrix and the untransformed 
non-colliding one, weighted with their respective probabilities, 

g' tot (At) = Prob(C At \g tot )SM (f?tot|CU t ) S* + Prob (C At \g tot )M(gtot\G At ) 

= sr 1 / 2 ^ 1 /^^ + etot - r 1/2 Qtotf 1/2 At . (5.i5i) 

Here, the complementary map M. y\Q A t) is fixed by the requirement that the 
state gtot should remain unchanged both if the collision probability vanishes, 
r = 0, and if the scattering has no effect, S = I. 

Focusing on the nontrivial part T of the two-particle S-matrix S = I + iT 
one finds that the unitarity of S implies that 

Im(T)=i(T-Tt)=iTtT. (5.152) 
Using this relation we can write the differential quotient as 
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Q ' tot{A 2 t gtot = Trvwva-rt - ^TtTr^^totr 1 / 2 (5.153) 
-rr^^r^jtT + 1 [t + Tt. rvv^r 1 / 



It is now easy to arrive at a closed differential equation. We trace out the en- 
vironment, assuming, in accordance with the Markov approximation, that 
the factorization g to t = p ® Pe is valid prior to each monitoring inter- 
val At. Taking the limit of continuous monitoring At — > 0, approximating 
Tr E ([Re(T), V 1 / 2 Q tot V' L / 2 ]) ~ Trg ([r 1 / 2 Re(T)r 1 / 2 , ^tot])) and adding the 
generator H of the free system evolution we arrive at [34] 



i [H,p] +iTr E ( r^Re^r 1 '-. f > rL 



/2 

Tr B (Tr 1 / 2 r 1 / 2 ^ 



2 

-^^([p^p^r^^tTr 1 / 2 ) . (5.154) 

This general monitoring master equation, entirely specified by the rate oper- 
ator T, the scattering operator S = I + iT, and the environmental state pe, is 
non-perturbative in the sense that the collisional interaction is nowhere as- 
sumed to be weak. It is manifestly markovian even before the environmental 
trace is carried out, and one finds, by doing the trace in the eigenbasis of 
Pe, that is has the general Lindblad structure (|5.79p of the generator of a 
quantum dynamical semigroup. The second term in (|5.154p . which involves a 
commutator, accounts for the renormalization of the system energies due to 
the coupling to the environment, just like (|5.148[) . while the last three lines 
describe the incoherent effect of the coupling to the environment. 

So far, the discussion was very general. To obtain concrete master equa- 
tions one has to specify system and environment, along with the operators 
T and S describing their interaction. In the following applications, we will 
assume the environment to be an ideal Maxwell gas, whose single particle 
state 

P gas = §ex P (5.155) 

is characterized by the inverse temperature /3, the normalization volume fl, 
and the thermal de Broglie wave length yl t h defined in (|5. 128(1 . 



5.4.3 Collisional Decoherence of a Brownian Particle 



As a first application of the monitoring approach, let us consider the "local- 
ization" of a mesoscopic particle by a gaseous environment. Specifically, we 
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will assume that the mass M of this Brownian particle is much greater than 
the mass m of the gas particles. In the limit m/M — > the energy exchange 
during an elastic collision vanishes, so that the mesoscopic particle will not 
thermalize in our description, but we expect that the off-diagonal elements 
of its position representation will get reduced, as discussed in Sect. 15.3.41 

This can be seen by considering the effect of the S-matrix in the limit 
m/M — > 0. In general, a collision keeps the center-of-mass invariant, and only 
the relative coordinates are affected. Writing So for the S-matrix in the center 
of mass frame and denoting the momentum eigenstates of the Brownian and 
the gas particle by |P) and \p), respectively, we have [35] 

S|P) \p) = I d 3 Q|P -Q)\p + Q) (^p ~^-P + Q|So|^p - 2£P) , 
J m M m M 

(5.156) 

where m» = Mm/ (M + m) is the reduced mass and Q is the transfered 
momentum (and thus the change of the relative momentum) . In the limit of 
a large Brownian mass we have m*/m — > 1 and m*/M — > 0, so that 

S\P)\p) -> J d 3 Q\P-Q)\p + Q){p + Q\S \p) [forM>m]. (5.157) 

It follows that a position eigenstate \X) of the Brownian particle remains 
unaffected by a collision, 

S|X)|VWe = \X) (V^^Soe^/*) |^„)e , (5.158) 

as can be seen by inserting identities in terms of the momentum eigenstates. 
Here, |V>in).E denotes an arbitrary single-particle wave packet state of a gas 
atom. The exponentials in (|5. 158[1 effect a translation of So from the origin 
to the position X, so that the scattered state of the gas particle |Vw )e 
depends on the location of the Brownian particle. 

Just like in Sect. 15.1.11 a single collision will thus reduce the spatial co- 
herences p (X, X') = (X\p\X') by the overlap of the gas states scattered at 
positions X and X' , 

p'(X,X') = p(X,X'){^\^) E . (5.159) 

The reduction factor will be the smaller in magnitude the better the scattered 
state of the gas particle can "resolve" between the positions X and X' . 

In order to obtain the dynamic equation we need to specify the rate oper- 
ator. Classically, the collision rate is determined by the product of the current 
density j — n gauS v Te i and the total cross section a (p re i), and therefore T should 
be expressed in terms of the corresponding operators. This is particularly sim- 
ple in the large mass limit M — > oo, where v IC \ — \p/m — P/M\ — > \p\ /m, so 
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that the current density and the cross section depend only on the momentum 
of the gas particle, leading to 

T = n gas — ff(p) ■ (5.160) 

m 

If the gas particle moves in a normalized wave packet heading towards the 
origin then the expectation value of this operator will indeed determine the 
collision probability. However, this expression depends only on the modulus of 
the velocity so that it will yield a finite collision probability even if the particle 
is heading away form the origin. Hence, for (|5.154[) to make sense either the 
S-matrix should be modified to keep such a non-colliding state unaffected, or 
T should contain in addition a projection to the subset of incoming states, 
see the discussion below. 

In momentum representation, p(P,P') = (P\p\P'}, equation (|5.154p as- 
sumes the general structural 

IP 2 - (P') 2 

+ J dP o dP' p(P o ,P' )M (P,P';P ,P' ) 
-\ J dP p(P ,P') J dP f M(P f ,P f ;P ,P) 
-\ J dP' p(P,P' ) J dP f M(P f ,P f ;P',P' ) (5.161) 
The dynamics is therefore characterized by a single complex function 



M (P, P'; P 0) P' ) = <P| tr gas (T^ 2 [|P ) (P' \ ® Pgas ] r^jt] \ P >) , 

(5.162) 

which has to be evaluated. Inserting the diagonal representation of the gas 
state (|5.155j) 

Pgas = ^—Jf- y"dpoAi(Po) \Po){Pa\ (5.163) 
it reads, with the choices (|5. 157ft and (|5.160[) for S and f, 



The second term in (|5.154p describes forward scattering and vanishes for mo- 
mentum diagonal pE- 
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M(P, P',P-Q,P'-Q') = J d Pl d PoA1 (p )<5 (Q+Pi -po)S (Q' + po - Pi) 



*3 



x^|Pok(Po)^l(Pi|T |Po)| 2 
m \l 

= S(Q-Q') /dpoM(Po) — |Pok(po) 

x^|(Po-Q|T |po)| 2 
=:^(Q-Q')M n (Q) . (5.164) 

This shows that, apart form the unitary motion, the dynamics is simply 
characterized by momentum exchanges described in terms of gain and loss 
terms, 

d tP (P,P') = ^ P2 ^P P(P,P')+ I dQp(P-Q,P'-Q)M in (Q) 

- P (P,P') J dQM- m (Q) . (5.165) 

We still have to evaluate the function M- ln (Q), which can be clearly inter- 
preted as the rate of collisions leading to a momentum gain Q of the Brownian 
particle, 

M in (Q) = ^ / d PoA1 (po) |po| o (po) ^-|(po - Q|T |po)| 2 . 

(5.166) 

It involves the momentum matrix element of the on-shell To-matrix, So = 
1 + 'iT , which, according to elastic scattering theory [35j . is proportional to 
the scattering amplitude /, 

The delta function ensures the conservation of energy during the collision. 
At first sight, this leads to an ill-defined expression since the matrix element 
(|5.167|) appears as a squared modulus in (|5.166p . so that the tree-dimensional 
integration is over a squared delta function. 

The appearance of this problem can be traced back to our disregard of 
the projection to the subset of incoming states in the definition (|5. 160|) of V. 
When evaluating Mj n we used the diagonal representation (|5.163[) for p gas in 
terms of (improper) momentum eigenstates, which comprise both incoming 
and outgoing characteristics if viewed as the limiting form of a wave packet. 
One way of implementing the missing projection to incoming states would be 
to use a different convex decomposition of pgas? 

which admits a separation 
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into incoming and outgoing contributions [36j . This way, M; n can indeed 
be calculated properly, albeit in a somewhat lengthy calculation. A shorter 
route to the same result sticks to the diagonal representation, but modifies the 
definition of S in a formal sense so that it keeps all outgoing state invariant 
The conservation of the probability current, which must still be guaranteed 
by any such modification, then implies a simple rule how to deal with the 
squared matrix element |36j , 

<^'«-^f (?-!)■ ( , 168) 

Here a{p) = J di?' \f(pn',pn)\ 2 is the total elastic cross section. With this 
replacement we obtain immediately 



Mn (Q) = ^J d Po A 4 (po) \f(po - Q, Po)| 2 S - 



Po (Po - Q) 



2 



(5.169) 



As one would expect, the rate of momentum changing collisions is determined 
by a thermal average over the differential cross section da/df2 = |/| 2 . 

Also for finite mass ratios m/M a master equation can be obtained this 
way, although the calculation is more complicated [37l [38]. The resulting 
linear quantum Boltzmann equation then describes on equal footing the de- 
coherence and dissipation effects of a gas on the quantum motion of a particle. 

The "localizing" effect of a gas on the Brownian particle can now be seen, 
after going into the interaction picture in order to remove the unitary part of 
the evolution, and by stating the master equation in position representation. 
From Eqs. (|5.165[) and (|5. 169[) one obtains 

d t p(X, X') = -F(X - X')p(X, X') (5.170) 

with localization rate [36j 

F(x) = j~ dv v(v) n gas v J (l - e'™^-^'^) 

x \f(mvn 2 , mvni)\ 2 . (5.171) 

Here, the unit vectors ni, n 2 are the directions of incoming and outgoing gas 
particles associated to the elements of solid angle dJ?i and dJ?2 and v (v) is 
the velocity distribution in the gas. Clearly, F (x) determines how fast the 
spatial coherences corresponding to the distance x decay. 

One angular integral in (|5.171j) can be performed in the case of isotropic 
scattering, f(p/,Pi) — f (cos (pf ,pi) ; E = p? / 2m) . In this case, 



17 In general, even a purely outgoing state gets transformed by S, since the definition 
of the S-matrix involves a backward time-evolution 1351. 
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Fig. 5.2. The localization rate (|5. 172ft describing the loss of wave-like behavior 
in a Brownian particle state saturates for large distances at the average collision 
rate. In contrast, the Caldeira-Leggett model predicts a quadratic increase beyond 
all bounds (dashed line), see (|5.126[l . This indicates that linear coupling models 
should be taken with care if time scales are involved that differ strongly from the 
dissipation time scale. 



2 



F(x) = J dvv(v) n gas v ^a(mv) — 2ir J d(cos0) / ^cos 8; E — -^v 2 

x sine ( 2 sin (§) ^M) } , (5.172) 

with sinc(a;) = sin(x)/a; and 9 the (polar) scattering angle. 

The argument of the sine function is equal to the momentum exchange 
during the collision times the distance in units of H. As \X — X'\ — > the 
sine approaches unity and the angular integral yields the total cross section 
a so that the localization rate vanishes, as required. At very small distances, 
a second order expansion in the distance x is permissible and one obtains a 
quadratic dependence [35] , such as predicted by the Caldeira-Leggett model, 
see Eq. (I5.127p . However, once the distance \X — X'\ is sufficiently large 
so that the scattered state can resolve whether the collision took place at 
position X or X' the sine function in (15. 172ft suppresses the integrand. It 
follows that in the limit of large distances the localization rate saturates, at 
a value given by the average collision rate F(oo) — (avn ga s}, see Figure I5T21 

Decoherence in this saturated regime of large separations has been ob- 
served, in good agreement with this theory, in molecular interference exper- 
iments in the presence of various gases [40]. The intermediate regime be- 
tween quadratic increase and saturation was also seen in such experiments 
on momentum-exchange mediated decoherence, by studying the influence of 
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the heat radiation emitted by fullerene molecules on the visibility of their 
interference pattern [IT] . 

As a conclusion of this section, we see that the scattering approach permits 
to incorporate realistic microscopic interactions transparently and without 
approximation in the interaction strength. The results show clearly that linear 
coupling models, which imply that decoherence rates grow above all bounds, 
have a limited range of validity. They cannot be judged by their success 
in describing dissipative phenomena. Frequent claims of "universality" in 
decoherence behavior, which are based on these linear coupling models, are 
therefore to be treated with care. 

5.4.4 Decoherence of a Quantum Dot 

As a second application of the monitoring approach, let us see how the dy- 
namics of an immobile object with discrete internal structure, such as an 
implementation of a quantum dot, gets affected by an environment of ideal 
gas particles. For simplicity, we take the gas again in the Maxwell state 
(|5.155|) , though different dispersion relations, e.g. in the case of phonon quasi- 
particles, could be easily incorporated. The interaction between system and 
gas will be described in terms of the in general inelastic scattering amplitudes 
determined by the interaction potential. 

In the language of scattering theory the energy eigenstates of the non- 
motional degrees of freedom are called channels. In our case of a structureless 
gas they form a discrete basis of the system Hilbert space. In the following, 
the notation la), not to be confused with the coherent states of Sect. I5.3T41 
will be used to indicate the system eigenstates of energy E a . In this channel 
basis, p a p = (a\p\(3), the equation of motion (|5.154| takes on the form of a 
general discrete master equation of Lindblad type, 

E a + E a - Ep — £p \ jixao/3o 

OtPafj = -Paf) + 2^ P* a Po M a 

-\ E P"o?> E M T - \ E p«po E M 77° ■ ( 5 - 173 ) 

"0 7 0o 7 

The real energy shifts e a given below describe the coherent modification of 
the system energies due to the presence of the environment. They are due to 
the second term in (|5.154p and are the analogue of the Lamb shift (|5.148|) 
encountered in the weak coupling calculation. The incoherent effect of the 
environment, on the other hand, is described by the set of complex rate 
coefficients 

Kf° = H Tr £ (TW 2 [|oo)03b| ® p gas ] rVa-rt) 1/3) . (5.174) 

In order to calculate these quantities we need again to specify the rate op- 
erator f. In the present case, it is naturally given in terms of the current 
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density operator j = n gas p/m of the impinging gas particles multiplied by 
the channel-specific total scattering cross sections a (p, a) , 



V = y \a)(a\<8> n gas — ^ (p, a) . (5.175) 



Like in Sect. 15.4.31 this operator should in principle contain a projection to 
the subset of incoming states of the gas particle. Again, this can be accounted 
for in two different ways in the calculation of the rates (|5.174|) . By using a 
non-diagonal decomposition of p gas , which permits to disregard the outgoing 
states, one obtains 18 ! 

Kf° = Xlf ^ I dp dpoA* (po) f aao (p,Po) 

xr m (p,p )s(^^+E a -E ao ) , (5.176) 



with the Kronecker-like factor 



a f3 J 1 ^ E a — E ao — Ep — Ep a . . 

Xa P ' \0 otherwise. 1 J 

The energy shifts are determined the real parts of the forward scattering 
amplitude, 

e a = -27rfi 2 ^ f dp of i(po)Re[f aa (p ,po)}. (5.178) 
m J 

Some details of this calculation can be found in [33] . Rather than repeating 
them here we note that the result (|5.176|) can be obtained directly by using the 
diagonal representation (|5. 163ft of p gas and the multichannel-generalization 
of the replacement rule (|5. 168ft , 

(27T/1) 3 f xTf° faao (P,Po)f*00o fogg) 

x^ p2 ^ P ° +E a -E a )j .(5.179) 

The expression for the complex rates simplifies further if the scattering ampli- 
tudes are rotationally invariant, f aao (cos(p,po) >E =p^/2m). In this case 
we have 



For the special case of factorizing interactions, Hint = A ® Be, and for times 
large compared to all system time scales this result can be obtained rigorously 
in a standard approach [42] , by means of the "low density limit" scaling method 
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T u f3 _ a Pt 



Kp = XZ J dw («) %asVout (v) 2n j ^ d (cos 6) 

(TTh \ / TTh \ 

cosB-E^-v 2 ) f* m (cos8;E=-v 2 ) (5.180) 

with v (w)the velocity distribution like in (|5.172[) . and 



«outW = \/f 2 --(£ a --E a „) (5.181) 
V m 

the velocity of a gas particle after a possibly inelastic collision. 

This shows that limiting cases of (|5.173p display the expected dynamics. 
For the populations p aa it reduces to a rate equation, where the cross sections 
c QQ0 (E) = 271-/ d(cos#) \ f aao (cos6;E)\^ for scattering from channel a to 
a determine the transition rates, 



oca 



dvu (v) n sas v ont (v)a aao (t^ 2 ) ■ (5.182) 



In the case of purely elastic scattering, on the other hand, i.e. for M^f* = 
MapSaotobpPoi the coherences are found to decay exponentially, 

dt\ Pa p\ =-7a| stlc |/M • (5.183) 

The corresponding pure dephasing rates are determined by the difference of 
the scattering amplitudes, 

.clastic 



7 elas txc = 7r j dw( V )n sas V out (v) j ^ d (COS 0) 

(TTh \ / TTh \ 

co S 0;-v 2 ) - f pf} (cos6;-v 2 ) . (5.184) 

As one expects in this case, the better the scattering environment can dis- 
tinguish between system states \a) and \/3) the more coherence is lost in this 
elastic process. 

In the general case, the decay of off-diagonal elements will be due to a com- 
bination of elastic and inelastic processes. Although little can be said without 
specifying the interaction, it is clear that the integral over \f aa — in 
(|5.184p , a "decoherence cross section" without classical interpretation, is not 
related to the inelastic cross sections characterizing the population transfer, 
and may be much larger. In this case, the resulting decoherence will be again 
much faster than the corresponding relaxation time scales. 



5.5 Robust States and the Pointer Basis 



We have seen that, even though the decoherence predictions of linear coupling 
models has to be taken with great care, the general observation remains 
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valid that the loss of coherence may occur on a time scale 7j cco that is 
much shorter the relaxation time 7 -1 . Let us therefore return to the general 
description of open systems in terms of a semigroup generator C, and ask 
what we can say about a general state after a time t which is still small 
compared to the relaxation time, but much larger than the decoherence time 
scale. From a classical point of view, which knows only about relaxation, the 
state has barely changed, but in the quantum description it may now be well 
approximated by a mixture determined by particular projectors Pi, 



This set of projectors {P^}, which depend at most weakly on t, is called 
pointer basis or set of robust states [H]. It is distinguished by the fact 
that a system prepared in such a state is hardly affected by the environment, 
while a superposition of two distinct pointer states decoheres so rapidly that 
it is never observed in practice. 

We encountered this behavior with the damped harmonic oscillator dis- 
cussed in Sect. 15.3.41 There the coherent oscillator states remained pure un- 
der Markovian dynamics, while superpositions between (macroscopically dis- 
tinct) coherent states decayed rapidly. Hence, in this case the coherent states 
P a = \a)(a\ can be said to form an (over-complete) set of robust states, 
leading to the mixture 



with appropriate measure /i. 

The name pointer basis is well-fitting because the existence of such robust 
states is a prerequisite for the description of an ideal measurement device in a 
quantum framework. A macroscopic — and therefore decohering — apparatus 
implementing the measurement of an observable A is ideally constructed in 
such a way that macroscopically distinct positions of the "pointer" are ob- 
tained for the different eigenstates of A. Provided these pointer positions of 
the device are robust, the correct values are observed with certainty if the 
quantum system is in an eigenstate of the observable. Conversely, if the quan- 
tum system is not in an eigenstate of A, the apparatus will not end up in a 
superposition of pointer positions, but be found at a definite position, albeit 
probabilistically, with a probability given by the Born rule. 

The main question regarding pointer states is, given the environmental 
coupling or the generator £, what determines whether a state is robust or 
not, and how can we determine the set of pointer states without solving the 
master equation for all initial states. It is fair to say that this issue is not 
fully understood, except for very simple model environments, nor is it even 
clear how to quantify robustness. 




-1 



-1 



(5.185) 



t 




(5.186) 



274 Klaus Hornberger 



An obvious ansatz, due to Zurek [SI US], is to sort all pure states in the 
Hilbert space according to their (linear) entropy production rate, or rate of 
loss of purity, 

a t S n „[p] - -2tr(p£(p)) . (5.187) 

It has been called "predictability sieve" since the least entropy producing and 
therefore most predictable states are candidate pointer states [BJ. 

In the following, a related approach will be described, following the pre- 
sentation in [3, 46J. It is based on a time evolution equation for robust states. 
Since such an equation must distinguish particular states from their linear 
superpositions, it is necessarily nonlinear. 



5.5.1 Nonlinear Equation for Robust States 

We seek a nonlinear time evolution equation for robust pure states Pf which, 
on the one hand, preserves their purity, and on the other, keeps them as close 
as possible to the evolved state following the master equation. 

A simple nonlinear equation keeping a pure state pure is given by the 
following extension of the Heisenberg form for the infinitesimal time step, 

Pt+st = P t + 5t Q[A t , P t ] + [P f , [P t , B t ]]) , (5.188) 

where A and B are hermitian operators. In fact, the unitary part can be 
absorbed into the nonlinear part by introducing the hermitian operator 
Xf = — i[Af, Pt] + Bt. It "generates" the infinitesimal time translation of the 
projectors (and may be a function of Pt), 

P t+ st = P t +5t[P u [P u X t ]\ . (5.189) 

With this choice one confirms easily that the evolved operator has indeed the 
properties of a projector, to leading order in St, 

P\+st = Pt+st (5.190) 

and 

(Pt+stf = Pt+st +0(St 2 ) . (5.191) 
The corresponding differential equation reads 

dtPt = ?t+St ~ Pt = [Pt, [Pt, X t ]] . (5.192) 

To determine the operator X t one minimizes the distance between the time 
derivatives of the truly evolved state and the projector. If we visualize the 
pure states as lying on the boundary of the convex set of mixed states, then 
a pure state will in general dive into the interior under the time evolution 
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generated by C. The minimization chooses the operator X t in such a way that 
P t sticks to the boundary, while remaining as close as possible to the truly 
evolved state. 

The (Hilbert-Schmidt) distance between the time derivatives can be cal- 
culated as 

||£(Pt)-d t p t ||! s =tr[(z-[p t ,[p t ,x t ]]; 

=z 

= tr(z 2 -2(Z 2 P t -(ZP t ) 2 )) 

+2tr((Z-X) 2 P t -((Z-X)P t n . (5.193) 



We note that the first term is independent of X, whereas the second one is 
non- negative. With the obvious solution Xt = Z = C(Pt) one gets a nonlinear 
evolution equation for robust states Pt, which is trace and purity preserving 

BE 

d t Pt= [P t! [P t ,£(P t )]] . (5.194) 

It is useful to write down the equation in terms of the vectors |£) which 
correspond to the pure state P t = 

« e io = [Aio«i)-mio«i)io]io. (5-195) 

V « ' 

"decay rate" 

If we take C to be of the Lindblad form (|5 . T9[) the equation reads 



« = ^H|0 + £tJ< L £>*(U-<L 



; J-^( L ^ L_<L ^ U; 



\0 



~{H)e\S) . (5.196) 
in 

Its last term is usually disregarded because it gives rise only to an additional 
phase if (H)^ is constant. The meaning of the nonlinear equation (|5.196| is 
best studied in terms of concrete examples. 



5.5.2 Applications 

Damped Harmonic Oscillator 

Let us start with the damped harmonic oscillator discussed in Sect. 15 .3741 By 
setting H = huia^a and L = a (|5. 196[) turns into 

d t \0 - -Wa|0 + 7 (V ) € (a - <a>*) - \ (a f a)?)) 10 • (5.197) 
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Note that the hrst term of the non- unitary part vanishes if |£) is a coherent 
state, i.e. an eigenstate of a. This suggests the ansatz |£) = \a) which leads 
to 



d t \a) 



7, ,2 
— a 
2 1 1 . 



(5.198) 



It is easy to convince oneself that this equation is solved by 



= |a e- iw *- 7 * /2 > = e-l a 'l/V' a ' 



10} 



(5.199) 



with at — ao exp (— iuit — "ft/ 2). It shows that the predicted robust states are 
indeed given by the slowly decaying coherent states encountered in Sect. 15. 3.31 



Quantum Brownian Motion 

A second example is given by the Brownian motion of a quantum particle. 
The choice 

H = — and L = (5.200) 
2m A th y ' 

yields a master equation of the form (|5.117j) but without the dissipation term. 
Inserting these operators into (|5. 196|) leads to 

m = 2^ 10 " 7 ^ [(x_ <x> « )2 " ((x ~ < x >e) 2 >«M> ■ (5.201) 

The action of the non-unitary term is apparent in the position representation, 
£(2) = (x At positions x which are distant from mean position (x)^ as 

compared to the dispersion cr^(x) = ((x— (x)^) 2 )^ 2 the term is negative 
and the value £(x) gets suppressed. Conversely, the part of the wave function 
close to the mean position gets enhanced, 

(x\£) = l su PP ressed if \ x ~ (*k I > CT « W (5 202) 

v lw [enhanced if \x - (x) f | < <r ? (x) . v ' ; 

This localizing effect is countered by the first term in (|5.201|) which causes 
the dispersive broadening of the wave function. Since both effects compete 
we expect stationary, soliton-like solutions of the equation. 

Indeed, a Gaussian ansatz for |£) with ballistic motion, i.e. (p)^ = po, 
(x)^ = x +pot/m 7 and a fixed width cr^(x) = <r solves (|5.201[) provided [44] 

*-±^*-{*£5fT • < 5 ' 203 > 

see (|5.128[) . As an example, let us consider a dust particle with a mass of 10 fig 
in the interstellar medium interacting only with the microwave background of 
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T = 2.7 K. Even if we take a very small relaxation rate of 7 = 1/(13.7 x 10 9 y), 
corresponding to the inverse age of the universe, the width of the solitonic 
wave packet describing the center of mass is as small as 2 pm. This sub- 
atomic value demonstrates again the remarkable efficiency of the decoherence 
mechanism to induce classical behavior in the quantum state of macroscopic 
objects. 
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